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SUmiAEX 


The  Dredged-Naterial  Disposal  Management  Model,  D2M2,  is  a  simulation- 
optimization  model  designed  to  provide  information  required  to  answer  the 
following  disposal  management  questions: 

1.  If  a  specified  long-term  operation  policy  is  followed  for  an 
existing  or  proposed  dredged-material  disposal  system,  what  is  the 
final  status  of  the  system,  given  the  initial  conditions,  system 
physical  and  economic  characteristics,  and  dredged-material  volume 
estimates?  What  is  the  cost  of  following  the  operation  policy? 

Must  additional  disposal-site  capacity  be  provided? 

2.  What  is  the  least-costly  long-term  operation  policy  for  an  existing 
or  proposed  disposal  system,  given  initial  conditions,  system 
physical  and  economic  characteristics,  and  dredged-material  volume 
estimates? 

3.  What  is  the  least-costly  method  to  provide  the  additional  volume 
required?  Should  expired  leases  be  extended  or  should  new  sites  be 
acquired? 

4.  If  new  sites  are  to  be  acquired,  what  is  the  least-costly 
combination  of  sites?  What  is  the  least-costly  sequence  for 
acquiring  these  sites? 

5.  Are  alternative  material-management  alternatives,  such  as 
transferring  material  or  resting  sites,  cost  effective? 

The  information  is  provided  by  formulating  and  solving  a  mathematical 
programming  model  in  which  the  disposal  system  is  represented  as  a  network, 
as  illustrated  by  Fig.  1.  The  dredging  sites  and  disposal  sites  are 
represented  by  nodes  of  the  network.  These  nodes  are  connected  by 
capacitated  arcs  which  represent  the  facilities  for  transportation  of 
material  in  the  system.  A  unit  cost  is  assigned  to  "flow”  each  arc, 
representing  the  average  cost  of  moving  material  within  the  system.  For 
analysis  of  multiple-period  operation,  multiple  networks  are  developed  and 
linked  in  time.  Information  for  answering  the  capacity  expansion  questions 
is  provided  by  systematically  formulating  and  solving  network-flow 
programming  problems  which  represent  alterative  expansion  schemes;  the 
least-costly  scheme  is  identified  as  the  best  scheme.  Details  of  model 
formulation  are  presented  in  Appendices  I  and  II. 

Input  required  to  use  the  program  includes  (1)  an  estimate  of  quantity 
of  material  dredged  at  each  site  each  period,  (2)  a  description  of  physical 
and  economic  characteristics  of  each  disposal  site,  (3)  a  description  of 
physical  and  economic  characteristics  of  each  material-transportation 
facility,  and  (4)  specification  of  any  required  movement  of  material  within 
the  system.  A  detailed  description  of  input  is  presented  as  Appendix  IV. 


1 


a)  Disposal  system 


b)  Arc-node  representation  of  disposal  system 


FIG.  1.  -  Network  Representation  of  Disposal  System 


BACK6K0UND 


The  Corps  of  Engineers  has  been  responsible  for  maintenance  of  the 
navigable  waterways  of  the  United  States  since  1824.  The  maintenance 
includes  excavation  and  disposal  of  the  sediment  deposited  in  the  waterways. 
Current  common  practice  is  to  excavate  the  material  with  a  mechanical  or 
hydraulic  dredge  and  to  transport  it  to  a  disposal  site  either  by  pumping 
through  a  pipeline  or  by  carrying  the  material  to  the  site  in  barges  or  in 
hoppers  on  the  dredge.  The  disposal  site  may  be  an  offshore  site  selected  to 
minimize  interference  with  navigation,  or  the  disposal  site  may  be  a 
contained  upland  site.  Contained  disposal  sites  are  natural  or  man-made 
ponding  areas  into  which  the  dredged  material  is  pumped  or  lifted.  In  the 
disposal  site,  water  gradually  drains  and  evaporates  from  the  dredged 
material,  and  the  solids  densify  and  consolidate.  The  rate  of  dewatering, 
densifying,  and  consolidating  can  be  increased  by  surface  trenching,  wicking, 
surcharging,  or  pumping.  Detailed  descriptions  of  these  techniques  and  other 
technical  aspects  of  dredged-material  management  are  presented  in  Reference  1 
and  in  associated  reports  of  the  Corps*  Dredge  Material  Research  Program. 

Management  of  the  long-term  operation  of  a  dredged-material  disposal 
system  requires  selection  of  excavating  and  transporting  equipment, 
allocation  of  disposal-site  capacity,  selection  of  appropriate  disposal-site 
management  practices,  and  identification  of  capacity  expansion  schemes.  Due 
to  the  complexity  of  the  long-term  problem,  equipment  selection  and  capacity 
allocation  generally  is  addressed  only  for  the  short-term,  with  equipment  and 
sites  selected  for  minimum  cost  at  the  time  the  dredging  is  performed. 
Physical  and  environmental  limitations  may  constrain  this  selection  and 
allocation.  Likewise,  problems  of  disposal-site  operation  and  of  capacity 
expansion  generally  are  addressed  with  heuristic  rules  as  problems  arise. 

PROGRAM  CAPABILITIES 

The  Dredged-material  Disposal  Management  Model  (D2MZ)  is  a 
simulation-optimization  model  for  systematic  analysis  of  long-term  operation 
and  expansion  of  a  disposal  system.  Application  of  D2M2  requires  estimates 
of  volumes  dredged,  descriptions  of  the  existing  and  potential  disposal 
sites,  and  description  of  the  dredging  and  transporting  facilities.  With  the 
model,  system  disposal  capacity  expansion  alternatives  can  be  analyzed,  and 
the  minimum-cost  disposal-site  acquisition  and  lease  extension  schedule  can 
be  determined.  This  is  accomplished  by  evaluating  automatically  the  present 
value  of  alternative  sequences  of  acquiring  new  sites  and  extending  leases  to 
identify  the  least-costly  expansion  policy.  The  cost  of  any  alternative 
capacity  expansion  plan  is  considered  to  be  the  sum  of  site  acquisition  cost, 
fixed  operation,  maintenance  and  repair  cost,  lease  cost,  and  operation 
cost.  To  determine  the  operation  cost,  D2M2  includes  the  capability  to 
identify  the  minimum-net-cost  short-term  operation  policy  for  any  specified 
system.  This  is  accomplished  by  formulating  and  solving  a  mathematical 
programming  problem  that  represents  the  problem  of  allocating  efficiently  the 
available  capacity.  If  desired,  this  portion  of  the  program  may  be  used 
without  the  capacity  expansion  evaluation  portion. 
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Disposal-sit*  consolidation  rates,  containment  dike  heights,  and  other 
characteristics  of  existing  and  proposed  disposal  system  components  are 
specified  by  the  model  user.  Thus,  management  schemes  that  involve  changes 
in  these  parameters  may  be  evaluated  by  systematic  variation  and  re-execution 
of  the  model. 

CmPUTATIONAL  TBCHMIQUIS 

Program  D2M2  employs  network-flow  programming  and  a  heuristic 
enumeration  scheme  to  evaluate  disposal-system  management  alternatives. 
Network-flow  programming  is  used  for  evaluation  of  the  operation  of  a 
specified  disposal  system.  En’meration  is  used  for  selection  of  the 
least-costly  capacity-expansion  scheme. 

Network-flow  Programming.  -  Program  D2M2  models  the  characteristics  of 
a  system  of  dredge  sites,  disposal  sites,  and  material  transportation 
facilities  as  a  network.  This  network  includes  nodes  that  represent  the 
available  disposal  sites  and  dredge  sites.  These  nodes  are  connected  by  arcs 
that  represent  transportation  linkages,  material-transfer  facilities, 
material-reuse  capabilities,  and  material  storage  in  disposal  sites.  The 
flow  of  material  through  these  arcs  represents  the  transporting,  or  storing 
material  within  the  disposal  system.  A  detailed  description  of  this  network 
formulation  is  presented  in  Reference  2.  A  reprint  of  this  reference  is 
included  as  Appendix  I. 

A  unit  cost  is  associated  with  the  flow  in  each  arc;  the  objective  of 
the  solution  algorithm  is  to  determine  the  allocation  of  flow  to  the  various 
network  arcs  to  minimize  the  sum  of  the  product  of  flow  in  each  arc  and  the 
corresponding  ''ost.  The  unit  costs  assigned  to  the  network  arcs  are  the 
discounted  'on.t  costs  of  storing,  transporting,  or  transferring  material  or 
the  negative  of  the  unit  benefit  of  reusing  material.  A  specialized 
network-flow  programming  algorithm  is  used  to  determine  the  minimum-cost 
flows.  This  algorithm  is  described  by  Jensen  and  Barnes  (3). 

Annotated  examples  of  applications  of  the  network  model  are  presented  in 
Appendix  III. 

Enuauration.  -  The  least-costly  scheme  for  acquisition  of 
user-specified  capacity  expansion  options  is  determined  in  D2M2  by 
enumeration  of  a  limited  number  of  the  possible  schedules.  Determination  of 
the  cost  of  each  schedule  is  accomplished  by  computing  the  sum  of  the  present 
value  of  the  acquisition  costs  and  lease-renegotiation  costs  and  the  present 
value  of  operation  cost  of  the  disposal  system  with  the  expansion  sites 
available.  This  operation  cost  is  determined  with  the  previously  described 
network  model  of  the  disposal  system. 

The  heuristic  enumeration  procedure  incorporated  in  D2H2  provides  for  a 
well-structured,  systematic  search  of  the  site  acquisition  and  lease 
renegotiation  options.  Enumeration  begins  with  evaluation  of  the  total  cost 
if  all  capacity  expansion  sites  are  acquired  in  the  earliest  period  allowed 
and  if  all  leases  are  renegotiated.  This  capacity  expansion  scheme  is 
adjusted,  based  on  analysis  of  disposal  site  utilization  during  the  period  of 
analysis.  After  each  adjustment,  the  network  model  of  the  disposal  system  is 
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altered  accordingly,  and  the  least-costly  operation  policy  is  found.  The 
cost  of  acquisition  and  lease  renegotiation  for  the  new  expansion  scheme  is 
determined  and  is  added  to  the  operation  cost.  The  heuristic  rules  are  again 
applied,  and  the  process  is  repeated.  This  process  is  described  in  detail  in 
Appendix  II.  With  careful  application  of  the  procedure,  acceptable 
alternative  capacity  expansion  schemes  can  be  identified  with  reasonable 
computational  effort.  One  such  application  is  presented  in  Appendix  III  as 
example  five. 

DATA  SBQUIBSnniTS 

Physical  and  economic  characteristics  of  the  dredged-material  disposal 
system  must  be  defined  by  the  user.  A  detailed  description  of  the  input  is 
included  as  Appendix  IV. 

Bach  disposal  site  must  be  identified,  and  the  location  must  be  defined 
as  shown  in  Pig.  2.  The  distances  between  disposal  and  dredge  sites  are 
determined  from  the  specified  locations  and  are  used  to  select  the 
appropriate  unit  transportation  cost.  The  initial  conditions  of  the  site, 
the  maximum  allowable  rate  of  addition,  and  an  average  wet-to-dry  ratio  for 
the  site  must  be  specified.  Constraints  on  drying  the  site  must  be  defined, 
and  any  fixed  disposal  rates  must  be  specified.  The  storage-elevation- 
surface  area  relationship  must  be  defined  to  allow  complete  reporting  of 
simulated  disposal-site  operation.  If  capacity-expansion  costs  are  to  be 
computed,  the  cost  of  acquiring,  leasing,  and  maintaining  the  site  must  be 
specified. 

Bach  dredging  site  must  be  identified,  and  the  location  must  be  defined 
consistent  with  definition  of  disposal-site  location.  The  estimated  volume 
of  material  removed  each  period  must  be  specified,  and  sites  must  be 
identified  in  which  the  material  can  be  disposed.  The  facilities  for 
dredging  and  transporting  material  to  the  disposal  sites  must  be  defined. 

The  dredging  and  transporting  facilities  available  within  the  system 
must  be  described.  A  unit-cost  vs.  distance  function  must  be  specified  for 
each  transportation  facility;  the  cost  of  moving  material  from  dredging  sites 
to  disposal  sites  is  determined  with  this  function  and  distances  computed 
from  user-specified  locations. 

HABDUAKE  AMD  SOFTWASE  REQUISBinDITS 

Program  D2M2  is  written  in  ANSI  standard  FORTRAN  IV.  The  program  was 
developed  originally  on  a  CYBER  175  computer  but  has  been  implemented  and  is 
maintained  by  HEC  staff  on  a  Harris  500  computer.  Execution  of  D2M2  requires 
use  of  a  random  access  input-output  package;  such  a  package  is  available  on 
most  computers.  Compilation  of  D2H2  with  the  Harris  computer  requires  four 
central- processor  seconds.  Execution  of  the  example  problems  of  Appendix  II 
requires  57,  193,  182,  238,  and  716  seconds,  respectively,  and  1.2  million 
(octal)  24-bit  words  of  virtual  memory. 
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Prosraa  D2M2  uses  fiv*  scrttch  file*.  Thtse  fllai  «r«  aialgned  as 
follows : 

Flla  Idantif ication  Usa 


TAPES 

User  input  file 

TAPE6 

Output  file 

TAPBll 

File  for  storage  of  capacity-expansion 

Iteration  results. 
Input  file. 

Also  used  for  restart 

TAPE12 

Intermediate  file 
copied. 

to  which  user  Input  file  Is 

TAPE99 

Direct-access  file 

for  network  storage 

Allocation  of  available  computer  memory  to  the  arrays  of  program  D2M2  Is 
problem  dependant  and  Is  performed  automatically  by  the  program.  When  the 
dimensions  of  the  disposal  system  and  the  period  of  analysis  are  defined,  an 
Internal  algorithm  subdivides  the  siemory  as  necesssary  to  store 
characteristics  of  the  disposal  sites,  dredge  sites,  and  transportation 
links.  Thus  operation  of  a  system  with  many  disposal  sites  and  dredge  sites 
can  be  analyzed  for  a  few  time  periods,  or  operation  of  a  small  system  can  be 
analyzed  for  many  time  periods.  The  limitations  on  disposal  system  size  and 
period  of  analysis  are  presented  In  Appendix  IV  for  the  distributed  Harris 
500  version  of  the  progam;  users  stay  modify  program  dimensions  to  Implement 
the  program  on  other  computers . 

IBBOB  IDBMTIPICATKMi 

As  the  Input  for  program  D2M2  Is  read,  the  user-specified  values  are 
compared  with  expected  ranges  or  allowable  values.  If  an  error  Is  detected, 
a  message  is  printed.  These  messages,  which  are  self-explanatory,  are  listed 
In  Table  1.  Errors  are  classified  as  fatal  or  non-fatal.  If  a  fatal  error 
Is  discovered,  program  execution  terminates.  If  a  non-fatal  error  Is 
detected,  a  default  or  maximum  value  is  used  instead  of  the  user-specified 
value,  a  message  Indicating  this  Is  printed,  and  execution  continues. 

ORIGIN  OP  PEOGBAH 

Program  D2N2  was  developed  originally  In  the  Hydrologic  Engineering 
Center  (HEC)  with  financial  and  technical  assistance  from  the  Philadelphia 
District,  Corps  of  Engineers.  Subsequent  modifications  were  funded  by  the 
Philadelphia  District  and  the  Dredging  Divison  of  the  Water  Resources  Support 
Center  through  the  Waterways  Experiment  Station,  Corps  of  Engineers.  The 
software  was  developed  by  Rochelle  Barkin  and  the  network  model  and  the 
branch-and-bound  algorithm  were  formulated  by  David  Ford  and  Darryl  Davis, 
Chief,  Planning  Analysis  Branch.  Dr.  Ford  managed  program  development  and 
supervised  preparation  of  this  manual.  The  subroutines  for  solution  of  the 
network-flow  programming  problem  were  provided  by  Quentin  Martin  of  the  Texas 
Department  of  Water  Resources;  Dr.  Martin  provided  assistance  with 
application  of  these  subroutines.  Brian  Heverin  and  William  W.  L.  Lee  of  the 
Philadelphia  District  assisted  with  ntodel  testing  and  provided  valuable 
technical  guidance.  Bill  S.  Elchert  was  Director  of  the  HEC  during  model 
development . 


7 


Table  1.  -  Error  Messages 


Fatal  Error  Messages 


Message 


User  Action 


TOO  MANY  TRANSPORTATION  LINKS,  DREDGE  SITES. 
DISPOSAL  SITES,  OR  TIME  PERIODS.  <Equation 
for  program  maxlmum>.  CALCULATED  =  <value 
calculated  for  user  defined  system> 

MAXIMUM  =  <program  maximum> 


Reduce  number  of  transportation 
links,  dredge  sites,  disposal 
sites,  or  time  periods.  See 
input  description  for  maximum. 


TOO  MANY  NODES  IN  MODEL.  <Equation  for  Reduce  number  of  disposal  sites, 
program  maximum>.  CALCULATED  =  <user  value>,  dredge  sites  or  time  periods,  as 
MAXIMUM  =  <program  maximum>  indicated  by  equation. 


<id>  CARD  EXPErtED.  <id>  CARD  READ. 


Consult  input  guide  for  sequence 
of  cards. 


LATEST  ACQUISITION  DATE  OMITTED  FOR 
EXPANSION  SITE  <site  id>. 


CARD  NO.  <number>  IHE  VALUE  OF  lEXP 
(<user-specif ied  value>)  IS  INVALID 


Specify  latest  acquisition  date, 
field  3  of  SX  card. 


Replace  value  in  field  1  of 
specified  card  with  0,  1,  2,  or 
3. 


TOO  MANY  ARCS  ENCOUNTERED  IN  PERIOD 
<period  number> 


Reduce  number  of  dredge  sites, 
number  of  disposal  sites,  or 
number  of  transportation  links 
that  originate  or  terminate  at 
disposal  sites. 


INSUFFICIENT  AMORTIZATION  PERIODS 


CARD  NO.  <number>  HAS  INVALID  DREDGE  SITE 
ID  =  <alphanuroeric  identifier> 


The  latest  period  specified  for 
possible  acquisition  of  this 
capacity  expansion  site  will  not 
permit  proper  amortization. 
Decrease  number  of  periods 
required  (field  9  of  J1  card)  or 
increase  latest  period  (field  3 
of  SX  card) . 

Replace  identifier  on  specified 
card  with  one  of  identifiers 
defined  on  DI  cards. 
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MegsaKe 


User  Action 


CARO  NO.  <number>  HAS  INVALID  DISPOSAL 
SITE  ID  3  <alphanuinerlc  identifier> 


Replace  identifier  with  one 
defined  on  SI  cards. 


CARD  NO.  <number>  HAS  INVALID  TRANSPORTATION 
SITE  ID  >  <alphanumeric  identifier> 


Replace  identifier  with  one 
defined  on  TI  cards. 


THE  DISTANCES  OF  THE  DISTANCE-UNIT  COST 
FUNCTION  (<transportation-type  identifier>) 
ARE  NOT  INCREASING 


Re-arrange  values  on  TD  cards . 
Values  on  TC  card  must 
correspond. 


THE  UNIT  COSTS  OF  THE  DISTANCE-UNIT  COST 
FUNCTION (<transportat ion-type  identif ier>) 
ARE  NOT  INCREASING 


Re-arrange  values  on  TC  card. 
Values  on  TD  card  must 
correspond. 


Non-Fatal  Error  Messages 


THE  VALUE  OF  NPER  (<user-specif led  value>)  Reduce  NPER  (field  1  of  J1  card) 

ON  THE  J1  CARO  IS  GREATER  THAN  THE  MAXIMUM  or  accept  maximum. 

(<maximum  allowable  value>)  NPER  SET  EQUAL 
TO  MAXIMUM 


THE  VALUE  OF  DRATE  (<user-specif led  value>)  Accept  default  or  specify 
ON  THE  J1  CARD  IS  INVALID.  DRATE  SET  0<  DRATE  A  100. 

EQUAL  TO  DEFAULT  (1) 


THE  VALUE  OF  VMULT  (<uger-specif led  value>) 
ON  THE  DL  CARO  NO.  <card  number>  IS  INVALID. 
VMULT  SET  EQUAL  TO  DEFAULT  (1.0) 


Accept  default  or  specify 
VMULT  >  0. 


THE  VALUE  OF  WDRAT  (<user-8pecif led  value>) 
ON  THE  SL  CARO  NO.  <card  number>  IS  INVALID. 
WDRAT  SET  EQUAL  TO  DEFAULT  (1.0) 


Accept  default  or  specify 
wet-to-dry  ratio  (WDRAT)  >  0. 


THE  VALUE  OF  NREHND  (<user-8pecif led  value>)  Accept  default  or  specify 

ON  THE  SR  CARD  NO.  <card  number >  IS  INVALID.  NREHND  i  2. 

NREHND  SET  EQUAL  TO  DEFAULT  (2) 


REAOIN  —  NEXT  IDENTIFICATION  CODE  IS  Check  card  for  keypunch  error. 

INVALID,  RECORD  IGNORED 
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Dredged-Material  Disposal 
Management  Model 

By  David  T.  Ford,*  M.  ASCE 


Absthact;  To  identify  efficient  dredged-mateiul  disposal  management  strat¬ 
egies  for  the  Delaware  River  ruvigation  system  near  Philadelphia,  the  system 
operation  problem  is  formulated  and  solved  as  a  generalized  miiumum  cost 
network  flow  programming  problem.  This  formulabon  represents  material  sources 
and  available  disposal  sites  as  nodes  of  the  network  and  transportation  iiitks 
and  carry-over  storages  as  arcs.  The  dewatering,  consoUdation,  and  densifi- 
cation  of  dredged  material  is  modeled  with  an  arc  gain  factor,  thereby  allowing 
reduction  of  the  total  volume  of  material  within  the  network  but  requiring  use 
of  a  network-with-gains  algorithm  for  solution  of  the  operation  problem.  Ap¬ 
plication  of  the  model  defines  cost-efficient  dyiuunic  schemes  for  allocation  of 
material  to  available  disposal  sites.  A  generalized  computer  program  was  de¬ 
veloped  to  define  automatically  the  nodes,  arcs,  and  parameters  of  the  arcs  of 
the  network,  even  a  description  of  the  dredged-material  disposal  system. 
Structured  analysis  and  structured  programming  techiuques  were  used,  thus 
providing  a  clear  definition  of  the  computations  required,  the  order  in  which 
they  must  be  accomplished,  and  the  flow  of  data.  ‘Dus  software  development 
technique  reduces  the  effort  required  for  subsequent  modification  of  the  pro¬ 
gram  to  analyze  the  system  capacity-expansion  problem. 


Delaware  River  Disposal  Management  Problem 

Background. — ^The  Corps  of  Engineers  has  been  responsible  for  main¬ 
tenance  of  the  navigable  waterways  of  the  United  States  since  1824.  The 
maintenance  includes  excavation  and  disposal  of  the  sediment  deposited 
in  the  waterways.  Current  common  practice  is  to  excavate  the  material 
with  a  mechanical  or  hydraulic  dredge  (10)  and  to  transport  it  to  a  dis¬ 
posal  site  either  by  pumping  through  a  pipeline  or  by  carrying  the  ma¬ 
terial  to  the  site  in  barges  or  in  hoppers  on  the  dredge.  The  disposal  site 
may  be  an  offshore  site  selected  to  minimize  interference  with  naviga¬ 
tion  or  the  disposal  site  may  be  a  contained  upland  site.  Contained  dis¬ 
posal  sites  are  natural  or  manmade  ponding  areas  into  which  the  dredged 
material  is  pumped  or  lifted.  In  the  disposal  site,  water  gradually  drains  ' 
and  evaporates  from  the  dredged  material,  and  the  solids  densify  and 
consolidate.  The  rate  of  dewatering,  densifying,  and  consolidating  can 
be  increased  by  surface  trenching,  wicking,  surcharging,  and  pumping 
with  well  pants.  Detailed  descriptions  of  these  techniques  and  other 
technical  aspects  of  dredged-material  management  are  presented  in  Ref. 

3  and  in  .associated  reports  of  the  Corps'  Dredged  Material  Research 
Program. 

Management  of  the  long-term  operation  of  a  dredged-material  dis¬ 
posal  system  requires  selection  of  the  equipment  to  be  used  for  exca¬ 
vating  and  transporting  the  material  from  the  channel  to  the  disposal 
sites,  allocation  of  the  capacity  of  the  available  disposal  sites  to  satisfy 

‘Hydr.  Engr.,  The  Hydrologic  Engig.  Center,  U.S.  Army  Corps  of  Engr.,  Davis, 
Calif  95616. 

Note. — Discussion  open  until  June  1, 1984.  To  extend  the  dosing  date  one  month, 
a  written  request  must  be  filed  with  the  ASCE  Manager  of  Tech^cal  and  Profes¬ 
sional  Publications.  The  manuscript  for  this  paper  was  submitted  for  review  and 
possible  publication  on  October  14,  1982.  This  paper  is  part  of  the  Journal  of 
Water  Resources  Planning  and  Management,  Vol.  110,  No.  1,  January,  1984. 
©ASCE,  ISSN  0733-9496/84/0001-(X)51/$01.00.  Paper  No.  18526. 
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the  demand  for  storage  imposed  by  the  dredging  operation,  selection  of 
appropriate  disposal-site  management  practices,  and  identification  of  ca¬ 
pacity  expansion  schemes  if  the  system  capacity  is  exhausted  at  some 
time.  Due  to  the  complexity  of  the  long-term  problem,  equipment  se¬ 
lection  and  capacity  allocation  generally  is  addressed  only  for  the  short¬ 
term,  with  equipment  and  sites  selected  for  minimum  cost  at  the  time 
the  dredging  is  performed.  Physical  and  environmental  limitations  may 
constrain  this  selection  and  allocation.  Likewise,  problems  of  disposal- 
site  operation  and  of  capacity  expansion  generally  are  addressed  with 
heuristic  rules  as  problems  arise. 

Delaware  River  System. — The  Delaware  River,  Delaware  Bay,  and  as¬ 
sociated  tributaries  are  maintained  in  a  navigable  condition  by  the  Phil¬ 
adelphia  District,  Corps  of  Engineers.  Within  this  area,  shown  in  Fig.  1, 
23  Federal  navigation  projects  yield  approximately  8,100,000  cu  yd 
(6,200,000  m*)  of  dredged  material  annually.  Non-Federal  maintenance 
dredging  contributes  an  additional  3,400,000  cu  yd  (2,600,000  m^).  The 
material  is  disposed  in  21  containment  sites.  According  to  estimates  pub¬ 
lished  in  a  1979  study,  by  1999  all  these  sites  will  be  filled  or  unavailable 
due  to  lease  expiration  with  continued  maintenance  dredging  at  current 


RG.  1.— Delaware  River  Syetem 


rates  and  with  no  change  in  management  practices  (4).  This  in  turn  would 
mean  reduction  or  cessation  of  dredging  and  consequent  reduction  or 
cessation  of  navigation.  The  1979  study  identifies  a  number  of  manage¬ 
ment  alternatives  that  may  be  employed,  including: 

1.  Capacity  exp>ansion  alternatives:  (a)  Acquisition  of  new  upland  sites; 
(b)  open-water  disposal  of  dredged  material;  and  (c)  extension  of  leases 
on  sites. 

2.  Operation  alternatives:  (a)  IDewatering  of  disposal  sites;  (b)  increase 
in  containment  dike  height;  (c)  reuse  of  dredged  material;  (d)  reduction 
of  maintenance  dredging;  (e)  use  of  deposition  basins  to  reduce  shoal¬ 
ing;  (f)  reduction  of  sediment  erosion;  and  (g)  improvements  in  site 
management. 

MANAOeMENT  MOOEL  DESCRIPTION 

Model  Objective. — ^The  dredged-material  disposal  management  model 
was  developed  for  systematic  evaluation  of  and  comparison  of  alterna¬ 
tive  management  schemes.  With  the  model,  capacity  expansion  alter¬ 
natives  can  be  analyzed,  and  the  minimum-cost  combination  and  sched¬ 
ule  can  be  determined  for  new  site  acquisition  and  lease  extension.  Also, 
the  minimum-net-cost  operation  policy  for  any  specified  system  can  be 
determined.  This  policy  is  required  both  for  long-term  system  operation 
planning  and  for  solution  of  the  expansion  problem;  the  total  cost  of 
any  alternative  capacity  expansion  scheme  is  a  function  of  site  acquisi¬ 
tion,  lease  extension,  and  operation  costs.  The  minimum  operation  cost 
and  the  associated  operation  policy  are  determined  by  formulating  a 
mathematical  programming  model  that  represents  the  problem  of  allo¬ 
cating  efficiently  the  available  capacity.  Disposal-site  dewatering  rates, 
containment  dike  heights,  and  other  characteristics  of  the  disposal  sys¬ 
tem  are  specified  by  the  model  user,  so  management  schemes  that  in¬ 
volve  changes  in  these  parameters  are  evaluated  by  systematic  variation 
and  re-execution  of  the  model. 

Initial  development  of  the  disposal  management  model  is  linuted  to 
formulation  of  the  mathematical  programming  model  presented  herein 
for  analysis  of  operation  of  a  defined  system.  Ultimately  the  model  will 
be  expanded  to  address  the  capacity  expansion  problem,  using  a  branch- 
and-^und  algorithm  which  iteratively  enumerates  a  limited  number  of 
alternative  site  acquisition  and  lease  renegotiation  schemes,  evaluates 
the  efficiency  and  feasibility  with  the  operation  model,  and  identifies 
efficient  schemes  for  expanding  the  system.  The  branch-and-bound  pro¬ 
cedure  provides  rules  for  eliminating  from  consideration  many  costly  or 
infeasible  schemes  without  actual  evaluation  with  the  operation  model. 

Mathematical  Programming  Formulation. — The  mathematical  pro¬ 
gramming  formulation  of  the  dredged-material  disposal  system  opera¬ 
tion  problem  includes  continuity  constraints  for  material  sources  and  for 
disposal  sites,  transportation  link  and  disposal-site  capacity  constraints, 
and  carry-over  storage  constraints.  The  continuity  and  capacity  con¬ 
straints  define  the  operation  problem  for  each  period.  The  carry-over 
storage  constraints  relate  conditions  within  each  period,  yielding  a  multi¬ 
period  operation  problem.  Unit  costs  are  associated  with  transportation 
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and  disposal  of  dredged  material.  The  objective  is  to  minimize  die  total 
discounted  cost  of  system  operation.  This  formulation  is  similar  to  the 
solid-waste  disposal  model  formulated  by  Marks  and  Uebman  (13)  and 
to  the  wastewater  disposal  model  formulated  by  Brill  and  Nakamura  (2). 

A  continuity  constraint  is  included  for  each  material  source  and  for 
each  disposal  site  for  each  period  of  analysis.  The  form  of  the  equation 
for  each  material  source,  /,  for  each  peri^  T  is 

NDiSP 

2  FdJ.T)  -  V(/,T) . (1) 

/-I 

in  which  /  ==  index  of  disposal  sites;  NDISP  =  total  niunber  of  disposal 
sites;  F(IJ,T)  =  volume  of  material  transported  from  source  I  to  site  / 
in  period  T;  and  V{1,T)  -  total  volume  of  material  dredged  at  source  I 
during  period  T.  The  form  of  the  equation  for  each  disposal  site  /  for 
each  period  T  is 

[NDRG  -I  rNIXSP 

2  F(i,;,r)J  +  |  2  RT(/',/,T) 

/'»*/ 

-RT(/J',T)j  -RU(/,T)  =  S(/,r) . (2) 

in  which  NDRG  =  total  number  of  dredged-material  sources;  S(/,  T  - 
1)  =  volume  of  nuterial  stored  at  site  /  at  beginning  of  period  T  and  at 
the  end  of  period  T  -  1;  S(/,  T)  =  voliune  of  matei^  stored  at  site  /  at 
end  of  period  T;  RT  (/,/',  T)  *  volume  of  material  transferred  to  site  / 
from  site  /';  RT  (/',/,!)  *  volume  of  material  transferred  from  site  /  to 
site  }'  in  period  T;  RU  iJ,T)  *  volume  of  material  from  site  /  removed 
and  sold  for  reuse;  and  VF  (/)  =  an  average  volume-reduction  factor. 
The  volume  reduction  factor  reflects;  (1)  The  wet-to-dry  voliune  ratio  of 
the  dredged  material;  and  (2)  the  efficiency  of  the  disposal  site  manage¬ 
ment  practices.  The  wet-to^ry  volume  ratio  defines  the  average  volume 
of  dry  material  per  time  peri^  that  must  be  stored  at  the  disposal  site 
as  a  fraction  of  the  total  volume  of  material  in  situ.  E.g.,  a  wet-to-dry 
ratio  of  2.0  indicates  that  the  dredged  material,  when  wet,  will  occupy 
twice  the  volume  occupied  by  the  dried  material.  In  this  formulation  ffie 
volume  reduction  is  assumed  to  occur  within  one  period.  The  efficiency 
of  the  disposal  site  in  terms  of  achievement  of  the  reduction  depends 
on  the  site  management  techniques.  If  the  techniques  enq>lqyed  are  1(X)% 
efficient,  VF  (/)  will  equal  tna  reciprocal  of  the  wet-to-dry  ratio;  other¬ 
wise,  VF  (/)  1^  equal  the  product  of  this  reciprocal  and  the  estimated 
efficiency  of  the  dewatering  techniques  used  at  the  site.  Typical  values 
of  VF  (/)  range  from  0.50  to  1.00. 

The  total  volume  of  material  to  be  transported  to  or  from  a  disposal 
site  is  constrained  by  the  characteristics  of  ffie  pipeline,  hopper,  or  other 
device  used  for  transportation.  Likewise,  the  volume  of  material  depos¬ 
ited  at  a  site  each  period  T  is  conslrained  by  the  size  of  the  site.  These 
limitations  are  expressed  matitematicaUy  as 

F(I,/,T)sFMAX(J,7) . (3) 
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RTa/',T)^RTMAX(/J') . (4) 

RU(/,T)sRUMAX(/) . (5) 

S(/,T)  sSMAX(/) . (6) 


in  which  FMAX(/,/)  =  capacity  of  the  transportation  link  between 
dredged-material  source  I  and  disposal  site  /;  RTMAX (/,/')  =  capacity 
of  the  facilities  for  removing  material  from  disposal  site  /  and  transfer¬ 
ring  it  to  site  /';  RUMAXf/)  =  capacity  of  the  facilities  for  removing 
material  from  disposal  site  /  for  reuse;  and  SMAX(/)  =  storage  capacity 
of  disposal  site  /. 

In  addition  to  the  restrictions  on  transportation,  disposal-site  man¬ 
agement  practices  may  pose  a  limitation  on  the  rate  of  addition  of  “wet" 
material  to  the  site.  This  limitation  is  imposed  each  period  by  the  fol¬ 
lowing  constraint: 

NDRG 

2  F(/,/,T)  <ADDMAX(/) . (7) 

/-I 

in  which  ADDMAX(/)  =  maximum  allowable  volume  addition  per  period. 

The  operation  problem  is  to  determine  the  “best"  scheme  for  allocat¬ 
ing  the  material  dredged  each  period  to  the  available  sites  over  the  plan¬ 
ning  horizon.  The  efficiency  of  operation  is  defined  as  the  algebraic  sum 
of  the  present  value  of  costs  of  disposal  and  transportation  and  the  ben¬ 
efits  of  reuse.  Mathematically,  this  is  expressed  as 


NPERS  r  rNDRC  NDISP  'j 

z=  2  (i  +  R)-'  2  2  CF(/,/).F((,;,r) 

T=i  L I  i~i  ;-i  J 

'NDISP  NDRG  1  /-NDISP  NDISP  'v 

2  CS(/).  +1  2  2  CRT(/,;')*RT(;,/',r) 

.  /=!  /=I  J  I  /-I  /■•I  J 


{NDISI 

2 

/  =  ! 


CRU(/)*RU(/ 


-T)}] 


in  which  Z  =  the  present  value  of  system  net  benefits  for  the  period  of 
analysis;  R  =  discount  rate;  NPERS  =  number  of  time  intervals;  CF(/,/) 
=  unit  cost  of  transporting  material  from  dredge  site  I  to  disposal  site 
/;  CS(/)  =  unit  cost  of  adding  material  to  disposal  site  /;  CRT(/,/')  = 
unit  cost  of  removing  material  h'om  site  /,  transporting  to  and  disposing 
in  site  /';  and  CRU(/)  =  unit  benefit  of  reuse  of  material  from  site  /. 
These  costs  and  benefits  are  assumed  to  be  constant  over  time.  The  ob¬ 
jective  is  minimization  of  the  net  cost,  Z. 

Mathematical  Programming  Problem  Solution. — ^The  dredged-mate¬ 
rial  system  operation  model  as  presented  includes  linear  constraints  and 
a  linear  objective  function,  so  a  linear  programming  (LP)  algorithm  can 
be  used  to  determine  the  optimal  allocation  of  dredged  material  to  the 
available  sites.  However,  the  constraints  define  only  conservation  re¬ 
quirements  and  transportation  limitations,  and  the  costs  and  benefits  are 
functions  of  the  volume  of  material  transported  or  stored,  so  the  oper- 
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ation  problem  can  be  formulated  as  a  network-flow  programming  prob¬ 
lem.  In  a  network-flow  problem,  the  decisions  required  are  visualized 
as  flows  in  the  arcs  connecting  the  nodes  of  a  network,  and  the  objective 
is  to  choose  the  flow  in  each  arc  to  optimize  some  efficiency  measure, 
such  as  total  cost.  The  arcs  of  the  network  are  characterized  by  the  al¬ 
lowable  direction  of  flow,  the  maximum  and  minimum  amounts  of  flow 
that  can  pass  through  each  arc,  the  unit  cost  of  use  of  the  arc,  and  a 
gain  that  represents  the  fraction  of  flow  that  is  lost  (or  gained)  in  each 
arc.  The  constraints  are  limited  to  conservation  of  flow  at  the  nodes  of 
the  network  and  to  upper  and  lower  bounds  on  flows  in  the  arcs.  Al¬ 
gorithms  for  solution  of  the  network-flow  problems  are  more  efficient 
than  those  for  solution  of  the  general  LP  problems.  Ford  and  Fulkerson 
(8)  and  Jensen  and  Barnes  (12)  provide  detailed  descriptions  of  the  char¬ 
acteristics  of  network-flow  problems. 

The  network-flow  model  of  the  disposal  operation  problem  represents 
material  sources  and  available  disposal  sites  as  nodes  and  transportation 
links  and  carry-over  storage  as  arcs.  The  network  representation  of  a 
small  disposal  system  is  shown  as  Fig.  2.  Nodes  1  and  2  represent  the 
Pedricktown  North  disposal  site,  nodes  3  and  4  represent  the  Pedrick- 
town  South  site,  and  nodes  5  and  6  represent  Overflow  Site  1.  Nodes 
7,  8,  and  9  represent  the  Marcus  Hook,  Bellevue,  and  Cherry  Island 
dredge  sites,  respectively.  The  arcs  connecting  nodes  7,  8,  and  9  with 
nodes  1,  3,  and  5  represent  the  transportation  links  between  the  material 
sources  and  disposal  sites.  The  "flows"  in  these  arcs  represent  the  vol¬ 
umes  of  material  allocated  to  the  disposal  sites.  For  a  more  complex  sys¬ 
tem,  arcs  are  included  to  connect  each  source  with  each  site  available 
for  disposal  of  the  material  from  that  source.  An  upper  bound  is  im¬ 
posed  on  flow  in  these  arcs,  as  dictated  by  the  transportation  method  repre¬ 
sented. 

Removal  of  material  from  the  disposal  sites  for  reuse  is  represented 
by  an  arc  originating  at  the  disposal  site  node  and  terminating  at  a  node 


/ 


\ 


FIG.  2.— Single-Period  Example  of  Network  Representation 
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that  represents  the  point  of  sale  of  the  reused  material.  Material  re¬ 
moved  from  the  Pedricktown  South  site  and  transferred  to  the  Pedrick- 
town  North  site  is  represented  by  flow  in  the  arc  originating  at  node  4 
and  terminating  in  node  2. 

The  arc  originating  at  node  1  and  terminating  at  node  2,  the  arc  orig¬ 
inating  at  node  3  and  terminating  at  node  and  the  arc  originating  at 
node  5  and  terminating  at  node  6  are  included  as  a  computational  mech¬ 
anism  to  represent  the  drying  of  material  added  to  a  disposal  site  and 
to  limit  the  rate  of  addition  of  material  to  the  site.  The  gains  for  these 
arcs  are  the  volume  reduction  factors  [VF(/)  of  Eq.  2],  and  the  upper 
bounds  are  the  maximum  allowable  rates  of  disposal  [ADDMAXf/)  of 
Eq.  7]. 

Material  is  introduced  to  the  network  at  the  nodes  that  represent  the 
dredged-material  sites.  In  the  terminology  of  Jensen  and  Barnes  (12), 
these  volumes  are  node  external  flows;  the  quantity  of  flow  entering  the 
network  is  fixed. 

The  dashed  arcs  of  Fig.  2  represent  the  storage  of  material  in  the  sys¬ 
tem  disposal  sites.  The  flow  in  the  storage  arc  that  terminates  at  node 
2  represents  the  net  volume  of  dried  material  deposited  in  the  Pedrick¬ 
town  North  site  in  all  periods  prior  to  the  current  period  [S(/,  T  -  1)  of 
Eq.  2].  The  storage  arc  originating  at  node  2  represents  the  cumulative 
volume  of  dried  dredged-material  deposited  in  the  site  after  the  addition 
arikd  removal  of  material  in  the  period  shown  [S(/,  T  -  1)  of  Eq.  2].  When 
the  network  is  expanded  for  analysis  of  multiple-period  operation,  these 
storage  arcs  link  the  networks  that  represent  single-period  problems.  This 
is  shown  in  Fig.  3. 

A  unit  cost  is  associated  with  the  flow  in  each  arc;  the  objective  of  the 
solution  algorithm  is  to  determine  the  allocation  of  flow  to  the  various 
network  arcs  to  minimize  the  sum  of  the  product  of  flow  in  each  arc 
and  the  corresponding  cost.  The  unit  costs  assigned  to  the  network  arcs 
are  the  discounted  units  costs  of  storing,  transporting,  or  rehandling 
material  [CS(/),  CF (/,/),  and  CRT(/,/')  of  Eq.  8]  and  the  negative  of 
the  unit  benefit  of  reuse  [CRU(/)]. 

The  traditional  network  flow  programming  solution  algorithms,  such 
as  the  out-of-kilter  algorithm  (5,7,9)  were  developed  for  problems  in  which 


FIG.  3.— Multiple-Period  Network  Representation 
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all  gain  factors  are  unity  and,  thus,  are  not  applicable.  Consequently, 
for  solution  of  the  management  problem  as  formulated,  a  specialized 
network-with'gains  algorithm  is  employed.  This  algorithm  solves  the 
generalized  minimum-cost  network  flow  problem  with  any  nonnegative 
gain  factors  using  a  flow-augmentation  algorithm.  In  this  application, 
the  algorithm  begi.iS  with  flow  in  all  arcs  set  equal  to  zero.  The  mini¬ 
mum  cost  per  unit  of  additional  flow  to  each  node  of  the  network  and 
the  path  over  which  that  flow  may  be  obtained  is  determined.  The  total 
flow  through  the  network  is  increased  along  the  minimum-cost  path  tm- 
til  the  flow  in  one  or  more  arcs  in  the  path  exceeds  the  bounds.  This 
process  continues  iteratively  until  the  required  system  input  flows  are 
satisfied  or  a  maximum  possible  flow  through  the  network  is  obtained. 
This  algorithm  guarantees  achievement  of  a  feasible,  global  optimal  so¬ 
lution  if  such  a  solution  exist.  Additional  details  of  the  algorithm  are 
presented  by  Jensen  and  Bhaumik  (11)  and  by  Jensen  and  Barnes  (12); 
a  previous  application  of  the  algorithm  is  described  in  Ref.  14. 

Software  Oevelofment 

A  generalized  computer  program  was  developed  to  implement  the 
proposed  disposal-system  management  model  to  evaluate  alternative 
system  capacity  expansion  plans.  The  program  development  employed 
state-of-the-art  software  engineering  techniques,  including  structured 
analysis  and  structured  programming  (1). 

Structured  Analysis. — Structured  analysis  is  a  logical  process  for 
transforming  information  about  program  requirements  into  specifica¬ 
tions  for  the  program  that  is  to  be  developed.  This  approach  is  contrary 
to  usual  engineering  program  development  activities  in  which  everyone 
eagerly  gets  on  to  the  "real"  work — writing  code.  As  described  by  De¬ 
marco  (6),  the  structured  analysis  approach  has  the  following  charac¬ 
teristics:  (1)  It  yields  a  paper  model  of  the  program-to-be;  (2)  the  pro- 


RG.  4.— Top  Level  of  Control 
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gram  is  desigi\ed  in  a  top-down,  hierarchical  fashion  with  a  smooth 
progression  horn  abstract  definition  of  program  components  to  a  de¬ 
tailed  definition;  (3)  it  yields  a  set  of  connected  "mini-specifications"  of 
the  identified  program  components;  and  (4)  it  uses  diagrams  for  com¬ 
munication  of  ideas,  especially  between  the  program  user,  program  de¬ 
signer,  and  computer  system  analyst. 

Top-down  program  design  begins  with  the  establishment  of  firm  re¬ 
quirements  for  the  tasks  to  be  accomplished  by  the  program  and  with 
the  definition  of  data  required  to  accomplish  the  tasks.  The  overall  pro¬ 
gram  structure  (top-level)  is  then  defined,  with  progressive  refinement 
of  lower-level  components  of  the  program.  Fig.  4  shows  the  organization 
of  the  top  level  of  the  dredged-material  disposal  management  program. 


EVALUATE 
COST  OF  PLAN 


FIG.  6.— Coat  Evaluation  Componant 


FIQ.  7.— Op«ration*Co«l  Evaluation  Componant 


The  program  consists  of  an  "executive"  routine  controlling  an  "input," 
a  "process,"  and  an  "output"  routine.  Figs.  5  and  6  show  further  re¬ 
finement  of  the  process  component;  specification  of  the  other  compo¬ 
nents  is  refined  in  a  similar  manner.  Development  of  the  system  man¬ 
agement  model  was  planned  (and  funded)  for  completion  in  two  separate 
stages:  stage  1  includes  only  the  operation  model  development,  while 
stage  2  addresses  the  capacity  expansion  problem  in  more  detail,  hi  stage 
1,  several  of  the  components  shown  were  defined  only  in  concepti^ 
terms.  For  example,  detailed  specification  of  computational  techniques 
was  delayed  initially  in  the  case  of  component  3.1.  Nevertheless,  the 
data  transfers  and  the  required  resulto  of  execution  of  each  module  were 
defined. 

The  network-flow  programming  model  formulated  to  determine  the 
optimal  allocation  of  dredged  material  is  identified  as  component  3.3.4 
in  Figs.  6  and  7.  As  one  of  the  goals  of  program  development  is  to  pro¬ 
duce  a  management  model  usable  by  engineers  and  planners  who  are 
not  familiar  with  mathematical  programming  techniques,  component 
3.3.4. 1  is  included  here  to  translate  disposal-system  descriptive  data  into 
the  node-arc  representation.  The  resulting  generalized  minimum-cost 
network-flow  pr^lem  is  solved  with  code  included  in  component  3.3.4.2. 
Definition  of  this  component  is  further  refined  to  include  components 
of  the  network-with-gains  algorithm. 

Structured  Programming. — ^The  actual  computer  code  to  implement 
the  disposal  management  model  was  developed  firom  the  structured 
analysis  using  structured  programming  techniques.  Each  of  the  com¬ 
ponents  was  translated  into  one  or  more  subprograms  that  perform  in¬ 
dependently  single  tasks  required  for  solution  of  the  operation  problem. 
The  benefits  of  this  approach  are:  (1)  The  actual  development  time  is 
reduced  because  a  number  of  programmers  may  simultaneously  develop 
the  modules,  or  existing  code  may  be  used  easily;  (2)  complex  programs 
may  be  tested  in  parts,  with  each  module  verified  independently;  (3)  the 
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code  is  easier  to  understand  and  to  maintain;  (4)  the  resulting  code  is 
flexible  and  may  be  modified  by  changing  single  modules  indepen¬ 
dently;  and  (5)  documentation  of  the  code  is  easier.  Items  3,  4,  and  5 
are  significant  given  the  environment  within  which  the  computer  code 
described  here  will  be  used.  Although  a  single  application  motivated 
development,  application  to  other  disposal  operation  problems  is  likely 
and  despite  careful  program  design,  past  experience  indicates  a  frequent 
need  for  special-case  modification.  Often  these  modifications  must  be 
performed  by  someone  other  than  the  original  program  writer,  thus  the 
need  for  understandable  code. 

Application 

Operation  Evaluation. — The  operation  model  has  been  used  to  eval¬ 
uate  the  operation  of  various  existing  and  proposed  configurations  of 
the  Delaware  River  dredged-material  disposal  system,  including  the  sub¬ 
system  between  Philadelphia  and  the  sea.  As  modeled,  this  subsystem 
includes  19  dredge  sites  and  8  disposal  sites  (of  which  two  are  imaginary 
sites  for  overflow  if  the  system  capacity  is  insufficient).  Pertinent  data 
describing  the  disposal  sites  are  presented  in  Table  1.  Material  is  dredged 
at  average  annual  rates  shown  in  Table  2  and  is  transported  by  barge, 
hopper,  or  pipeline  to  the  disposal  sites.  The  dredging  and  transporting 
costs  depend  on  the  machinery  used  and  the  distance  which  the  material 
is  transported;  costs  vary  from  $1.62  to  $25.00/cu  yd  and  are  shown  in 
Table  3.  The  unit  costs  of  placing  the  material  in  the  system  disposal 
sites  vary  from  $0.00  to  $0.50/cu  yd,  as  indicated  in  Table  1. 

The  operation  of  the  Philadelphia-to-sea  subsystem  was  analyzed  for 
50  yr  using  25  consecutive  2-yr  interv'als.  The  resulting  network  con¬ 
sisted  of  8^  nodes  and  3,709  arcs.  Time  required  for  definition  of  the 
network  parameters  and  for  solution  of  the  minimum-cost  optimization 
program  on  a  commercial  Cyber  175  computer  system  used  by  the  Corps 
was  approximately  59  CP  sec.  The  minimum  present-value  net  cost  for 
system  operation  with  an  annual  discount  rate  of  7-5/8%  is  $273  x  10*. 

Figs.  8-11  are  reproductions  of  portions  of  the  computer  program  out¬ 
put.  Fig.  8  is  a  summary  of  the  dimensions  of  the  system  to  be  analyzed. 


TABLE  1.— Dltposai  Sit*  Information 


Disposal  site 
(1) 

Capacity  remaining, 
in  cubic  yards 
(2) 

Wet-to-dry 

ratio 

(3) 

Disposal  cost, 
in  dollars  per 
cubic  yard 
(4) 

Artificial  Island 

16,500,000 

1.50 

0.32 

Overflow  Site  1 

99,000,000 

1.00 

0 

Overflow  Site  2 

99,000,000 

1.00 

0 

National  Park 

7,100,000 

1.50 

0.50 

Killcohook 

36,900,000 

1.50 

0.11 

Penns  Neck 

16,000,000 

1.50 

0.25 

Pedricktown  North 

21,700,000 

1.50 

0.16 

Pedricktown  South 

21,700,000 

1.50 

0.17 

Note:  1  cu  yd  =  0.765  1  acre  =  0.405  ha. 
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TABLE  2.— Ofdgtng  Rate.  In  Cubic  Yards  par  Yaar 


Orsdge  sits 

(1) 

Volunis 

(2) 

Eddystone 

12,300 

Chester 

890 

Marcus  Hook 

1,850,300 

Bellevue 

49,900 

Cherry  Island 

180,900 

Deepwater  Point 

1,402,900 

Bulkhead  Bar 

28,700 

Newcastle 

1,269,400 

Reedy  Island 

28,000 

Baker 

10,800 

Liston 

220,800 

Miah  Maull 

41,000 

Brandywine 

1,500 

W.  Horseshoe 

25,800 

Mifflin 

67,400 

Billingsport 

5,600 

Tinicunt 

43,200 

Upper  Philadelphia  Harbor 

6,500 

Lower  Philadelphia  Harbor 

181,500 

Note:  1  cu  yd  =  0.765  m*. 

TABLE  3.-~Ortdglng  ami  Transporting  CoaL  In  Dollars  par  Cubic  Ysrd 


DispoMi  Sit* 


Dredge  Site 
(1) 

Artificial 

Island 

(2) 

Overflow 
Site  1 
(3) 

Overflow 

Site2 

(4) 

National 

Park 

(5) 

KillCO- 

hook 

(6) 

Penna 

Nock 

(7) 

PjJHl 

rfOnCK* 

town 

South 

(») 

Eddystone 

6.99 

25.00 

25.00 

17.89 

» 

• 

2.83 

3.06 

Chester 

6.57 

25.00 

25.00 

* 

A 

• 

2.42 

2.65 

Marcus  Hook 

6.12 

25.00 

25.00 

17.01 

• 

1.99 

2.22 

Bellevue 

4.91 

25.00 

25.00 

* 

• 

2.47 

2.09 

1.86 

Cherry  Island 

• 

25.00 

25.00 

15.40 

• 

2.08 

2.48 

2.25 

Deepwater  Point 

3.57 

25.00 

25.00 

.2.01 

1.95 

• 

3.23 

Bulkhead  Bar 

3.16 

25.00 

25.00 

1.61 

2.35 

• 

• 

Newcastle 

2.48 

25.00 

25.00 

• 

1.98 

3.02 

• 

A 

Reedy  Island 

1.62 

25.00 

25.00 

• 

284 

3.90 

* 

• 

Baker 

1.63 

25.00 

25.00 

• 

3.25 

4.29 

• 

• 

Liston 

2.22 

25.00 

25.00 

• 

3.85 

4.86 

• 

A 

Miah  Maull 

7.25 

19.25 

19.25 

• 

8.88 

9.94 

• 

Brandywine 

8.70 

15.24 

15.24 

• 

10.33 

11.36 

* 

• 

W.  Horseshoe 

8.56 

25.00 

25.00 

19.30 

* 

S 

4.37 

Mifflin 

8.09 

25.00 

25.00 

19.02 

S 

• 

3.93 

4.15 

Billingsport 

7.85 

25.00 

25.00 

18.78 

• 

• 

3.69 

3.93 

Tinicum 

7.23 

25.00 

25.00 

18.14 

S 

• 

3.07 

3.31 

Upper  Philadelphia 
Harbor 

9.69 

25.00 

25.00 

19.30 

• 

s 

5.50 

5.75 

Lower  Philadelphia 
Harbor 

9.24 

25.00 

25.00 

19.30 

• 

• 

5.29 

5.04 

Note:  *  Indicates  sites  not  linked  (1  cu  yd  *  0.765  m^). 
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Also,  tor  each  type  of  dredge,  a  function  relating  the  unit  excavating  and 
transporting  cost  to  distance  transported  is  presented.  Fig.  9  is  one  of  8 
disposal  site  reports.  In  this  report,  the  physical  and  economic  charac¬ 
teristics  of  the  disposal  site  are  summarized.  Site  acquisition  and  lease 
renegotiation  data  are  included  for  future  use  when  the  program  is  ex¬ 
panded  to  address  the  capacity  expansion  problem.  Fig.  10  is  one  of  19 
dredge  site  reports  for  this  system.  The  alternative  sites  for  disposal  of 
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material  are  shown,  and  the  capacity  and  the  unit  cost  of  excavating  and 
transporting  material  to  each  site  is  tabulated.  (The  unit  cost  is  deter¬ 
mined  from  the  appropriate  unit  cost  versus  distance  function.)  The  es¬ 
timated  volumes  of  material  to  be  removed  each  period  are  shown.  (Note 
that  the  value  shown  for  each  period  in  Fig.  10  is  twice  the  correspond¬ 
ing  value  from  Table  2  because  each  period  in  the  analysis  corresponds 
to  2  yr.) 

The  results  of  optimal  operation  of  the  dredged-material  disposal  sys¬ 
tem  are  presented  as  a  tabulation  of  material  added  to  each  disposal  site 
each  period  and  of  end-of-period  storage  and  the  corresponding  eleva¬ 
tion  and  surface  area.  Fig.  11  is  an  example  of  the  tabulation  for  one 
site.  The  results  of  the  minimum-cost  operation  are  summarized  in  Table 
4.  All  disposal  sites  except  the  overflow  sites  will  be  filled  by  the  end 
of  the  50-yr  model.  Overflow  Site  1  is  used  initially  in  the  ninth  2-yr 
interval,  indicating  that  system  capacity  falls  short  of  demand  within  18 
yr  of  the  first  year.  This  conclusion  is  significant  but  is  difficult  to  draw 
with  traditional  mass  analysis  techniques  because  of  the  complex  inter¬ 
connections. 

Systematic  Evaluation  of  Operation  Alternatives. — Evaluation  of  dis¬ 
posal  site  management  alternatives  is  accomplished  by  systematic  ap¬ 
plication  of  the  operation  model  with  variation  of  the  appropriate  input 
parameters.  For  example,  to  evaluate  the  cost  effectiveness  of  use  of 
trenching  devices  that  speed  the  drying  of  deposited  material  in  the  dis¬ 
posal  sites,  the  volume-reduction  factor,  VF(/),  and  the  maximum  al¬ 
lowable  volume  addition  per  period,  ADDMAX(/),  are  changed  to  re¬ 
flect  the  improvement  possible,  and  the  disposal  cost,  CS(J),  is  altered 
as  appropriate.  The  operation  problem  is  resolved  to  determine  the  least- 
costly  operation  scheme.  The  cost  of  the  trenching  machinery  is  added 
to  the  operation  cost  to  determine  the  total  system  cost.  This  total  cost 
is  compared  with  the  total  cost  without  the  trenching  devices;  if  the  cost 
is  less,  the  trenching  device  is  cost  effective. 

Any  management  alternative  can  be  evaluated  by  systematic  analysis 
with  the  operation  model,  if  the  improvements  attributable  to  that 
alternative  can  be  expressed  in  terms  of  the  volume-reduction  factor. 
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TABLf  4.— VoluiM  of  DrodgocMIolortel  CMopotod,  In  MHIion  CuMe  Yardo  botow 
Prying  


DOpOMiSto 


Dradg*  Site 

(1) 

Artificial 

Island 

(2) 

Ovarflow 
Slla  1 
(3) 

Ovarflow 

SNaZ 

W 

Nationai 

Park 

(5) 

KMoo- 

hook 

(6) 

Panna 

- «- 

(7) 

i- 

r^vOflGK* 

town 

Soulh 

(9) 

Eddystone 

0 

0.20 

0 

wrrm 

• 

• 

0.20 

0.17 

Chester 

0 

0.02 

0 

• 

* 

• 

0.01 

0.01 

Marcus  Hook 

0 

29.25 

0 

8.63 

• 

• 

27.25 

27.39 

Bellevue 

0 

0.90 

0 

• 

• 

1.00 

0 

0.60 

Cherry  Island 

• 

2.17 

0 

1.45 

• 

5.43 

0 

0 

Deepwater  Point 

0 

25.25 

0 

27.32 

17.58 

• 

0 

Bullhead  Bar 

0 

0.52 

0 

0.92 

0 

• 

• 

Newcastle 

14.92 

21.44 

0 

27.11 

0 

• 

• 

Reedy  Island 

0.95 

0.45 

0 

• 

0 

0 

• 

• 

Baker 

0.37 

0.17 

0 

0 

0 

• 

Liston 

7.51 

3.M 

0 

• 

0 

0 

• 

• 

MUh  MauU 

0.98 

1.07 

0 

• 

0 

0 

• 

• 

Brandywine 

0.05 

0 

• 

0 

0 

• 

• 

W.  Horseshoe 

0 

0.52 

0 

• 

0.36 

0.36 

Mifflin 

0 

1.35 

0 

0 

• 

» 

0.94 

1.08 

Billingsport 

0 

0.10 

0 

0.01 

• 

0.09 

0.08 

Tinicum 

0 

0.78 

0 

• 

• 

0.69 

0.60 

Upper  Philadelphia 
Haibor 

0 

0.13 

0 

0.01 

• 

0.10 

0.06 

Lower  Philadelphia 
Harbor 

0 

3.63 

0 

0.36 

2.90 

2.18 

Total 

24.75 

91.53 

0 

10.65 

55.35 

24.01 

32.54 

32.55 

Note:  *  Indicates  sites  not  linked  (1  cu  yd  »  0.765  m^). 


maximum  allowable  addition  per  period,  maximum  storage,  capacity  of 
transfer  facilities,  capacity  of  reuse  facilities,  or  the  time  at  which  facil¬ 
ities  are  available.  Techniques  that  alter  the  volume  of  time  distribution 
of  dredged  material  that  must  be  disposed  can  be  analyzed  in  the  same 
systematic  manner  because  the  volumes  are  specified  by  the  user  for 
each  period. 

Conclusions 

To  identify  efficient  dredged-material  disposal  management  strategies, 
the  system  operation  problem  can  be  formulated  and  solved  as  a  gen¬ 
eralized  minimum-cost  network  flow  programming  problem.  In  this  for¬ 
mulation,  the  material  sources  and  disposal  sites  are  represented  as  nodes, 
and  the  transportation  links  and  carry-over  storages  are  represented  as 
arcs.  The  network-with-gains  algoritlm  is  used  for  solution,  thereby  al¬ 
lowing  modeling  of  the  drying  of  material  in  the  disposal  sites. 

The  pro(>osed  disposal  system  operation  model  can  be  used  for  the 
evaluation  of  alternative  management  schemes  by  systematically  varying 
the  appropriate  model  parameters,  re-executing  the  model,  and  com¬ 
paring  the  net  operation  cost  to  determine  the  effectiveness  of  the  scheme. 
In  the  future,  this  model  will  be  linked  with  a  branch-and-bound  algo¬ 
rithm  to  identify  efficient  disposal  system  expansion  schemes. 
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A  generalized  computer  program  to  implement  the  proposed  opera¬ 
tions  model  was  developed  using  software  engineering  m^ods.  Struc¬ 
tured  analysis  techniques  were  used  to  define  the  program  require¬ 
ments,  and  structured  programming  was  used  to  transform  the  require¬ 
ments  into  executable  computer  code. 

The  dredge-material  disposal  management  model  has  been  used  suc¬ 
cessfully  to  evaluate  the  operation  of  the  Delaware  River  disposal  system. 
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ApnmoL  It— Notation 


The  following  symbols  are  used  in  this  paper: 

ADDMAX(/)  «*  maxiinum  allowable  volume  addition  per  period,  site 

Ji 

CF(J,/)  -  urut  cost  of  transporting  material  from  source  I  to 
disposal  site  J; 

unit  cost  of  rehandling  material  from  disposal  site  / 
to  site  /'; 

CRU(/)  »  unit  benefit  of  reuse  of  material  from  site /; 

CS(/)  «  unit  cost  of  adding  material  to  site  /; 

F(I,J,T)  s  volume  of  material  transported  from  source  /  to  site 
/  in  period  T; 

FMAX  (/,  / )  s  capacity  of  transportation  link  between  dredged-ma- 
terial  source  7  and  disposal  site  /; 

NDISP  »  number -of  disposal  sites; 

NDRG  »  number  of  dr^ged-material  sources; 

NPERS  *  number  of  the  periods; 

RT(/,/',T)  =  volume  of  material  transferred  from  site  /  to  site  7 
period  T; 

RTMAX(/,/')  =  capadtyoftransferfadlities  between  site/and  site/'; 
RU(/,  T)  =  volume  of  material  removed  from  site  /  and  sold  for 
reuse,  period  T; 

RUMAX(/)  «  capacity  of  reuse  facilities,  disposal  site  /; 

S(//T)  =  volume  stored  in  disposal  site /,  period  T; 

SMAX(J)  *  capacity  of  dispell  site  /; 

V(7,  T)  -  total  volume  of  material  dredged  at  source  7  during 
period  T; 

VF(/)  *  volume  reduction  factor,  disposal  site  /;  and 
Z  =  net  system  operating  cost. 

Subscripts 

7  »  dredged-material  source; 

7  »  disposal  site;  and 
T  =  time  period. 
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APPDIDIX  II 


D2112  CAPACITY  KZPANSION  ALGORITHM 


SUMMARY 

The  capacity  expansion  aiodule  of  program  D2M2  selects,  from  a 
user-specified  set  of  alternatives,  the  least-costly  plan  for  acquiring  or 
leasing  disposal  sites  to  expand  the  capacity  of  a  disposal  system.  The 
least-costly  plan  is  determined  by  systematic  enumeration  and  comparison  of  a 
limited  number  of  acquisition  and  lease-renegotiation  schedules.  The  cost  of 
each  is  the  present  value  of  acquisition  and  lease-renegotiation  cost  plus 
the  present  value  of  disposal-system  operation  cost.  The  operation  cost  is 
determined  by  formulating  a  network-flow  model  of  the  alternative  schedules 
and  solving  the  minimum-cost  flow  problem  to  determine  optimal  allocation  of 
the  available  space.  A  branch-and-bound  enumeration  procedure  is  used  to 
select  candidate  expansion  plans  for  evaluation. 

BRANCH-AND-BOUND-MBTHODS 

General  Properties.  -  Branch-and-bound  methods  are  intelligently 
structured  techniques  for  evaluating  and  comparing  solutions  to  optimization 
problems.  With  these  methods,  only  a  small  fraction  of  all  possible 
solutions  actually  are  enumerated.  The  remaining  solutions  are  eliminated 
from  consideration  through  establishment  of  progressively  stricter  bounds  on 
the  cost  or  benefit  of  solutions  (the  objective  function).  These  bounds  aid 
in  identification  of  inferior  solutions  without  requiring  extensive 
computation  to  evalute  those  solutions. 

Branch-and-bound  techniques  consist  of  two  basic  operations:  (1) 
dividing  the  set  of  possible  solutions  into  subsets  (branching),  and  (2) 
establishing  bounds  on  the  value  of  the  objective  function,  so  inferior 
solutions  can  be  eliminated  (bounding).  The  methods  involve  recursive 
application  of  these  operations,  with  specific  rules  for  the  operations 
depending  on  the  characteristics  of  the  problem. 

References.  -  The  general  characteristics  of  branch-and-bound  methods 
and  applications  of  the  methods  have  been  presented  in  the  management  science 
and  operations  research  literature.  Lawler  and  Wood  (6)  present  a  survey  of 
the  essential  features  of  branch-and-bound  methods  for  constrained 
optimization  problems  and  describe  application  to  integer  and  nonlinear 
programming  problems,  to  the  traveling-salesman  problem,  to  the  quadratic 
assignment  problem,  and  to  non-mathematical  programming  problems.  Mitten  (9) 
describes  the  general  properties  of  branch-and-bound  methods  and  presents,  in 
general  terms,  the  conditions  for  branching  and  for  bounding  the  results  of 
optimization  problems.  Garfinkel  and  Nemhauser  (5)  describe  branch-and-bound 
methods  applicable  to  integer  programming  problems. 

Branch-and-bound  methods  have  been  applied  in  resource  planning  to 
problems  of  sizing,  selecting,  sequencing,  and  scheduling  projects.  Marks 
and  Liebman  (8)  suggest  using  a  branch-and-bound  procedure  for  locating 
solid-waste  management  facilities.  Brill  and  Nakamura  (3)  and  Nakamura  and 
Brill  (11)  propose  application  of  a  branch-and-bound  method  for  generating 
systematically  alternative  plans  for  regional  wastewater  treatment  systems 
and  for  evaluating  these  alternative  plans.  Application  of  branch-and-bound 


29 


••thods  for  selectloa  of  the  optiaal  coablnetioa  of  discrete  slternetlves  Is 
addressed  by  Efroynson  and  Ray  (4),  by  Morin  (10)  and  by  Ball,  Bialas,  and 
Loocks  (1). 

The  branch-and-bound  procedure  incorporated  in  prosraa  D2M2  is  an 
adaptation  of  the  procedure  sussested  by  Lesso,  Hissaelblau,  Jensen,  and 
Shanausham  (7)  and  the  procedure  suggested  by  Bickel  (2).  These  references 
provide  a  auitheawtically  vigorous  explanation  of  the  procedure. 

CAPACITY  RZPAIISION  ALGORITHM 

The  branch-and-bound  algorithm  of  program  D2M2  identifies  the  least- 
costly  capacity  expansion  plan  by  (1)  defining,  for  each  site,  a  time  window 
during  which  the  site  may  be  acquired,  (2)  partitioning  this  tisw  window 
repeatedly,  (3)  determining  a  lower  bound  on  the  cost  if  the  sites  are 
acquired  within  the  partitions,  (4)  analyxing  the  utilization  of  each  disposal 
site  within  the  defined  partitions,  (5)  estiauiting  the  cost  of  alternative 
partitions  for  each  site,  (6)  reducing  the  partition  for  the  site  which  will 
yield  the  greatest  reduction  in  overall  cost,  and  (6)  repeating  the  process 
until  the  partition  for  each  site  is  reduced  to  a  single  period  or  until  an 
iteration  limit  is  exceeded. 

Branching.  -  The  branch-and-bound  procedure  partitions  the  set  of  all 
possible  capacity  expansion  plans  into  progressively  smaller  subsets  and 
compares  the  cost  of  these  subsets.  The  partitioning,  or  bounding,  is 
accoaqilished  using  an  acquisition  time  window.  For  each  expansion  site,  the 
program  user  defines  the  earliest  and  latest  acceptable  acquisition  periods. 
The  many  alternative  acquisition  options  represented  by  this  acquisition  time 
window  are  evaluated  by  partitioning  repeatedly  the  initial  window  into  small 
candidate  segments  and  evaluating  the  cost  of  each.  Bach  of  these  segments 
represents  a  capacity-expansion  plan  with  site  acquisition  sometime  between 
two  periods  within  the  user-specified  tinte  window  The  partitioning  continues 
until  the  candidate  segments  are  reduced  to  a  single  period;  this  defines  the 
acquisition  period. 

Fig.  3  shows  how  the  branch-and-bound  algorithm  might  partition  a 
user-specified  acquisition  tisie  window  for  evaluation.  Site  acquisition  is 
possible  betkreen  2000  and  2050.  For  analysis  of  system  operation  with  the 
site  acquired  after  2000  and  before  2050,  a  partit.ion  is  made  yi»lding  two 
candidates:  acquisition  after  2010  and  in  or  before  2050.  Analysis  of  these 
options  leads  to  partitioning  of  the  latter  candidate.  (The  option  of 
acquiring  the  site  between  2000  and  2010  will  be  evaluated  later.)  The 
partition  yields  two  candidate  plans:  acquisition  between  2010  and  2020  or 
acquisition  between  2020  and  2050.  Following  analysis  of  system  operation, 
the  option  of  acquiring  the  site  between  2011  and  2020  is  partitioned.  The 
resulting  candidates  specify  acquisition  between  2011  and  2015  or  between  2016 
and  2020.  The  earlier  acquisition  promises  to  provide  greater  cost  reduction, 
so  it  is  analyzed,  and  an  additional  partition  is  made.  In  this  partition, 
the  site  may  be  acquired  in  2011  or  in  or  after  2012  and  in  or  before  2015. 

If  the  first  option  is  feasible  a,  solution  to  the  capacity  expansion  probloa 
has  been  found. 


30 


If  lavaral  capacity  expansion  sites  are  available,  the  capacity 
expansion  plans  are  partitioned  based  on  combinations  of  the  acquisition 
periods.  For  example,  the  program  user  may  specify  that  site  A  can  be 
acquired  between  2000  and  2050,  and  site  B  can  be  acquired  between  1990  and 
2010.  The  solution  can  be  partitioned  into  candidates  in  which  A  is  acquired 
between  2000  and  2025  with  B  acquired  between  1990  and  2010  and  candidates  in 
which  A  is  acquired  between  2026  and  2050  with  B  acquired  between  1990  and 
2010.  The  first  solution  can  be  partitioned  further  into  solutions  in  which 
A  is  acquired  between  2000  and  2025  with  B  acquired  between  1990  and  1998  and 
solutions  with  the  same  acquisition  window  for  A  but  with  B  acquired  between 
1999  and  2010.  This  partitioning  continues  until  a  single  acquisition  period 
is  found  for  each  site. 

Evaluating  Cost.  -  The  true  cost  of  each  capacity-expansion  plan  is 
the  present  value  of  acquistion  and  lease-renegotiation  cost  plus  the  present 
value  of  disposal-system  operation  cost.  If  the  acquisition  period  for  each 
site  is  known,  this  cost  can  be  computed  by  discounting  appropriately  the 
acquisition  cost  from  the  period  of  acquisition,  adding  the  discounted  OMR 
cost,  and  adding  the  discounted  operation  cost.  The  discounted  operation 
cost  is  computed  by  formulating  a  network  model  and  solving  to  define  the 
least-costly  operation.  During  the  search  for  the  optimal  capacity-expansion 
plan  with  partitions  of  the  acquisition  time  window,  the  period  of 
acquisition  for  any  candidate  capacity-expansion  plan  is  known  only  to  fall 
between  the  earliest  and  latest  periods  irtiich  define  the  partition. 
Consequently  the  acquisition  and  lease  renegotiation  cost  cannot  be  computed 
exactly.  For  each  disposal  site  it  is  estimated  as  a  function  of 
disposal-site  utilization  as  follows: 

1.  The  annual  equivalent  of  the  acquisition  cost  is  determined  using  a 
user-specified  amortization  period.  The  annual  OMR  cost  is  added. 

2.  The  resulting  cost  is  divided  by  the  capacity  of  the  site,  yielding 
a  unit  cost  approximation. 

3.  This  unit  cost  is  added  to  the  user-specified  cost  of  disposing 
material  in  the  site.  The  present  value  of  the  unit  cost  is 
determined  for  each  period  between  the  earliest  and  latest  defined 
by  the  partition.  The  resulting  unit  costs  are  assigned  to  the 
network  model  arcs  representing  storage  in  the  disposal  site. 

The  total  cost  computed  using  the  estimate  described  always  under¬ 
estimates  or  is  exactly  equal  to  the  true  cost.  Lesso,  Himmelblau,  Jensen, 
and  Shanmugham  (7)  prove  this.  This  fact  is  critical  for  eliminating 
inferior  solutions  without  explicitly  evaluating  the  solutions.  If  two 
alternative  expansion  plans  are  compared,  the  plan  with  the  greatest  cost  is 
inferior,  regardless  of  the  size  of  the  partition.  Even  if  another,  smaller 
partition  of  the  inferior  plan  is  defined,  the  cost  will  not  decrease.  Thus 
all  solutions  in  which  the  site  is  acquired  between  the  periods  defining  the 
partition  can  be  eliminated  from  further  consideration.  Any  effort  to 
increase  the  earliest  period  or  to  decrease  the  latest  period  is  wasted. 

Bounding.  -  The  branch-and-bound  procedure  estimates  the  cost  of  each 
candidate  solution  and  compares  the  cost  with  an  upper  bound  on  the  total 
cost.  This  upper  bound  is  the  value  of  the  best  solution  found  in  the 
iterative  procedure  and  is  updated  when  better  solutions  are  found.  If  the 
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value  of  any  candidate  solution  exceeds  the  upper  bound,  the  candidate  auy  be 
eliminated  from  further  consideration.  Furthermore,  because  the  cost  of  the 
candidate  represents  a  lower  bound  for  solutions  that  might  result  from 
additional  partitioning,  all  additional  partitions  can  be  eliminated  also. 

Bstiauiting  Cost  Seduction.  -  Selection  of  the  partition  of  the 
acquisition  time  window  for  further  investigation  is  accomplished  by 
estimating  possible  cost  reduction.  Each  capacity  expansion  site  is 
considered  in  turn,  and  the  cost  of  the  unused  capacity  is  estimated.  Later 
acquisition  of  a  site  should  decrease  this  cost,  but  it  also  renders  the 
project  unusable  in  earlier  periods,  so  additional  cost  is  incurred  then. 

The  algorithm  computes  the  sum  of  the  possible  cost  decrease  and  the  cost 
increase  due  to  modification  of  the  acquisition  period.  The  site  that 
promises  the  greatest  cost  decrease  is  selected  for  further  partitioning  of 
the  solution.  The  partition  is  BMde  in  the  period  which  yields  this  decrease. 

Reiterating.  -  When  a  site  is  selected  for  partitioning  and  a  period 
is  selected  in  which  the  partitioning  is  to  occur,  the  latest  or  earliest 
period  for  acquisition  of  the  site  is  redefined  to  equal  the  selected 
period.  Thus  a  new  partition  of  the  acquisition  time  window  is  defined  for 
that  site;  the  partition  is  smaller  than  the  originally-defined  window.  The 
costs  are  then  re-estimated  and  the  evaluation  is  repeated.  The  procedure 
terminates  when  the  earliest  and  latest  acquisition  periods  coincide  for  each 
project  or  when  a  user-specified  iteration  limit  is  reached. 
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APPKNOIX  III 


02M2  EXAMPLES 

Program  02H2  provides  the  capability  to  model  a  variety  of  disposal 
systems  and  to  simulate  a  variety  of  system  operation  criteria.  This 
appendix  includes  five  examples  trfaich  illustrate  program  capabilities.  Each 
example  includes  a  description  of  the  problem,  with  input  to  and  output  from 
program  D2M2.  Key  items  of  the  input  and  output  are  numbered  for  convenience 
of  explanation. 

Example  one  illustrates  application  of  D2M2  to  analyze  the  continuous 
operation  of  an  existing  disposal  system.  In  example  two,  the  restart 
capability  of  the  program  is  used  to  analyze  operation  for  tiro  saccessive 
25-year  planning  periods.  Example  three  is  a  modification  of  example  one, 
with  provisions  for  **resting**  the  Edgewood  site  for  two  periods  after  each 
three  periods  of  availability.  Example  four  illustrates  simulation 
specification  of  a  portion  of  the  operation  policy.  The  capacity  expansion 
capability  of  the  program  is  illustrated  by  example  five. 

EXAMPLE  ONE;  ANALYSIS  OF  CONTINUOUS  OPERATION 

This  example  illustrates  application  of  program  D2M2  to  determine  the 
status  of  and  operation  cost  of  a  dredged-material  disposal  system  if  the 
least-costly  operation  policy  is  followed  for  SO  periods,  each  of  which  is 
one  year  in  duration.  The  system  consists  of  two  dredging  sites  from  which 
material  may  be  transported  to  four  disposal  sites,  as  illustrated  by  Fig. 

4.  Analysis  begins  in  1981.  The  ocean  disposal  site  and  the  Williamson 
Harbor  and  Edgewood  sites  are  available  in  1981  and  remain  available 
throughout  the  SO-year  analysis  period.  The  lease  on  the  Edgewood  site  is  to 
be  renegotiated  in  1999  at  a  cost  of  $10,000.  The  Williastson  Harbor  South 
disposal  site  is  not  available  until  1987,  when  it  is  acquired  at  a  cost  of 
$1,725,000.  This  cost  is  amortized  with  7-5/8X  per  period  discount  rate. 

The  present  value  of  all  costs  is  determined  using  the  same  discount  rate. 

Two  dredging  and  transporting  techniques  are  employed:  a  16-incb  pipeline 
dredge  and  a  government-owned  hopper  dredge.  The  unit  cost  of  dredging  and 
transporting  the  material  is  a  function  of  the  distance  from  the  dredging 
site  to  the  disposal  site. 

Reproductions  of  selected  portions  of  the  D2M2  output  are  included 
herein.  Key  items  are  numbered  for  convenience  of  explanation. 

The  program  banner  page  (item  1)  identifies  the  version  of  the  program 
and  the  latest  revision  of  the  users  manual.  The  date  and  time  of  program 
execution  are  also  shown. 

Item  2  is  a  list  of  all  user-provided  input,  exactly  as  read  by  the 
program.  The  cards  are  numbered  for  subsequent  reference.  The  input  for 
this  example  includes  three  title  cards  (cards  1,  2,  and  3),  and  a  J1  card 
(card  4}  for  system  dimension  specification  and  job  control.  The  values  on 
the  J1  card  indicate  that  50  periods  will  be  analyzed,  beginning  in  1981 
(fields  1  and  2).  Each  period  in  this  case  is  one  year,  but  the  calculations 
performed  do  not  require  or  assume  this  duration.  Fields  4,  5,  and  6  of  the 
J1  card  define  the  number  of  disposal  sites,  nuisber  of  dredge  sites,  and 
number  of  dredged-material  transportation  methods.  The  value  in  field  9 
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Williamson  Harbor 
Dredging  Site 
(WIUfHA) 


Williamson  Harbor  B 
Dredging  Site 
(WILIIBO) 


Ocean  Disposal  Site 
(OCEAN  1) 


Williamson  Harbor 
Disposal  Site  (WILHAR) 


Edgewood  Disposal  Site 
(EDGEWD) 


Williamson  Harbor  South 
Disposal  Site  (VJILSUD) 


FIG.  4  .  -  Disposal  System  for  Example  One 


specifies  15  periods  for  amortization  of  site  acquisition  costs;  the  discount 
rate  of  7-5/8  percent  per  period  is  specified  in  field  3.  Identification  of 
the  amortization  period  is  required  to  allow  proper  computation  of  the 
present  value  of  acquisition  costs  for  comparison  of  alternative  schemes.  A 
set  of  cards  TI,  TD,  and  TC  cards  (cards  5-10)  is  included  for  each  of  the 
two  dredging  techniques.  These  cards  identify  the  dredging  technique  and 
describe  the  distance-cost  function  for  each.  The  values  of  unit  cost  as  a 
function  of  distance  are  used  with  computed  distance  to  estimate  appropriate 
unit  costs  between  interconnected  sites.  Linear  interpolation  is  used,  so 
the  values  specified  must  be  selected  accordingly.  Cards  11-16  define  the 
characteristics  of  the  Ocean  disposal  site,  cards  17-22  define  the 
characteristics  of  the  Williamson  Harbor  South  disposal  area,  cards  23-28 
describe  the  Williamson  Harbor  disposal  area,  and  cards  29-34  describe  the 
Edgewood  site.  Fields  1  and  2  of  the  SL  cards  define  the  location  of  each 
site;  the  program  uses  these  to  compute  distances.  The  off-river  distance 
for  the  ocean  site  (field  2  of  card  12)  is  specified  as  an  arbitrary,  large 
value  to  yield  a  high  cost,  thus  discouraging  ocean  disposal.  The 
acquisition  status  of  each  site  is  defined  by  the  value  in  field  1  of  the  SX 
cards  (cards  13,  19,  25,  and  31).  The  ocean,  Williamson  Harbor,  and  Edgewood 
sites  are  added  to  the  system  in  1981,  as  specified  by  the  value  in  field  2 

of  the  SX  cards.  Acquisition  costs  fcr  these  sites  are  ignored  in  this 

analysis,  so  these  are  omitted  from  the  SX  cards.  The  ocean,  and  Williamson 
Harbor  sites  are  available  throughout  the  analysis,  as  specified  in  field  4 

of  the  SX  cards.  In  the  case  of  the  Edgewood  site,  the  lease  is  to  be 

renegotiated  in  1999,  as  specified  by  the  value  in  field  4  of  card  31.  That 
site  is  then  available  throughout  the  remainder  of  the  analysis.  The 
Williamson  Harbor  South  site  is  added  to  the  system  in  1987  (field  2  of  card 
19)  and  is  available  throughout  the  analysis.  The  Williamson  Harbor  dredging 
site  is  described  by  DI,  DL,  and  DT  cards  (cards  35-39),  and  the  volume  of 
material  dredged  in  each  period  is  specified  on  DV  cards  (cards  40-45).  For 
this  example,  the  total  volume  is  specified  for  each  period,  so  six  DV  cards 
are  required  with  nine  values  specified  on  each.  Cards  46-56  describe  the 
Williamson  Harbor  B  site  and  define  the  volume  dredged  each  period. 

Each  input  error  detected  by  the  program  is  identified  (item  3).  The 
sequence  number  of  the  card  and  the  variable  which  is  specified  incorrectly 
are  shown.  If  a  default  value  is  used  rather  than  the  specified  value,  the 
default  value  is  shorn. 

The  basic  data  summary  (item  4)  shows  the  values  specified  by  the  user 
on  the  J1  card.  Item  5  is  the  summary  of  dredge  and  transportation  data. 

The  disposal  site  input  reports  summarize  the  data  for  each  disposal 
site  in  turn.  The  report  for  the  Williamson  Harbor  South  disposal  area 
comprises  items  6-17.  Similar  reports  are  printed  for  each  site  but  are  not 
shown  here.  Item  6  is  the  specified  location  of  the  disposal  site  in  terms 
of  main-stem  river  distance  and  off-main-stem  distance.  Any  origin  for 
measurement  may  be  used,  but  all  river  distances  must  be  specified  relative 
to  the  same  origin.  The  off-main-stem  distance  may  be  specified  as  a 
positive  or  negative  distance,  depending  on  the  location  to  the  right  or  left 
of  the  main  river.  The  off-river  transportation  type  also  is  shown.  If  this 
field  is  blank,  the  transportation  type  is  assumed  to  be  the  same  as  the 
transportation  type  used  on  the  main-stem.  Item  7  is  the  disposal  site 
storage  data,  including  the  initial  volume  of  material  in  the  containment 
site,  the  wet-to-dry  ratio  (volume  of  wet  material  divided  by  volume  of 
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•aterial  after  drying) ,  the  maxiauia  allowable  rate  of  addition  of  wet 
material,  and  the  unit  cost  for  addition  of  material  to  the  disposal  site. 

For  the  Williamson  Harbor  South  Disposal  area,  no  maximum  allowable  rate  is 
specified  (field  6  of  card  19),  so  the  site  capacity  is  used  as  the  default 
value.  The  site  acquisition  cost  (in  1981  dollars),  acquisition  date,  and 
the  expansion  indicator  (field  1  of  card  19)  are  shown  (item  8).  The  present 
value  of  the  acquisition  cost  is  determined  and  is  added  to  the  system 
operation  cost.  If  the  site  is  leased  and  a  lease-termination  date  is 
specified,  this  is  shown  (item  9).  Lease  renegotiation  costs  are  also  shown 
here.  The  present  value  of  the  renegotiation  costs  is  determined  and  is 
added  to  the  sum  of  acquisition  cost  and  operation  cost.  Data  on  reuse  of 
dried  material  from  the  disposal  site  are  summarized  (item  10),  and  transfer 
data  are  reported  if  provided  (item  11).  No  transfer  of  material  from  this 
site  is  possible,  so  no  data  are  shown  here.  If  transfer  is  possible,  the 
potential  ultimate  disposal  sites  and  the  transportation  techniques  are 
identified.  Finally  the  capacity-elevation-area  function  for  the  disposal 
site  is  shown  (item  12).  Linear  interpolation  is  used  with  these  values,  so 
they  should  be  specified  accordingly. 

The  dredge-site  input  report  for  the  Williamson  Harbor  dredging  site 
(items  13,  14,  and  IS)  shows  the  volumes  of  material  dredged,  by  period,  and 
the  techniques  available  to  accomplish  the  dredging  and  transporting  of 
material  to  the  disposal  sites.  A  similar  report  is  prepared  for  each  site, 
but  these  are  not  shotm  here.  The  dredging  site  location  (item  13)  is 
specified  in  terms  of  main-stem  river  distance  and  off-main-stem  distance, 
using  a  constant  origin  for  main-stem  locations  and  positive  or  negative 
distances  for  off-main-stem  distances.  The  material  volume  multiplier  is  a 
convenience  to  allow  scaling  of  all  specified  volumes  for  this  disposal  site 
prior  to  analysis  of  system  operation.  If  desired,  this  may  be  used  to 
account  for  the  difference  of  volume  in  situ  and  volume  to  be  transported. 

The  transportation  link  data  summary  (item  14)  identifies  disposal  sites  in 
which  material  from  this  dredging  site  may  be  disposed  and  the  techniques  for 
and  unit  cost  of  moving  the  material  to  the  sites.  Finally  the  material 
volumes  are  tabulated  (item  15);  these  are  the  user-supplied  values,  scaled 
by  the  volume  multiplier. 

The  volume  summary  report  (item  16)  shows  the  total  volume  by  period, 
before  drying  and  consolidating,  that  must  be  disposed  within  the  specified 
system.  The  dimensions  of  the  resulting  network-flow  programming  model  are 
shown  ( item  17) . 

The  capacity  expansion  iteration  log  serves  as  a  cost  summary  in  this 
application  of  D2H2.  Only  a  single  execution  on  the  model  is  required;  this 
is  shown  as  iteration  0  (item  18).  The  disposal  sites  are  identified  and 
pertinent  acquisition  and  lease  dates  are  shotm.  The  sum  of  present  value  of 
site  acquisition,  lease  renegotiation,  and  operation,  maintenance,  and 
replacement  (OMR)  costs  for  all  sites  is  computed  and  shown  (item  19).  The 
minimum-cost  present  value  of  system  operation  cost  is  shown,  as  determined 
by  solution  of  the  network-flow  programming  problem  (item  20).  This  cost 
represents  the  cost  of  material  transportation  and  disposal,  less  the  revenue 
from  reuse.  The  total  net  cost  (item  21)  is  the  total  of  acquisition,  lease 
renegotiation,  (HIR,  and  operation  costs  for  all  sites  for  the  period  of 
analysis . 
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The  results  of  optimal  operation  of  the  dredsed-material  disposal  system 
are  presented  as  reports  of  the  status  of  and  changes  to  each  disposal  site. 
These  reports  Include  a  tabulation  of  material  added  to  the  disposal  site, 
either  from  dredging  sites  or  from  other  disposal  sites,  a  tabulation  of 
material  removed  from  the  site  for  transfer  or  reuse,  and  a  tabulation  of  the 
containment  site  storage,  elevation,  and  surface  area.  An  example,  the 
report  for  the  Williamson  Harbor  South  site,  comprises  Items  22-25.  Similar 
reports  are  presented  for  all  sites,  but  are  not  included  here. 

The  Williamson  Harbor  South  disposal  site  may  be  used  for  containment  of 
material  from  the  Williamson  Harbor  dredging  site  and  from  the  Williamson 
Harbor  B  dredging  site.  The  volume  of  material  transported  from  each  is 
shown  and  the  discounted  transportation  cost  Is  tabulated  (Item  22).  The 
total  for  all  dredging  sites  is  shorn  (item  23),  and  the  total  for  all 
periods  is  computed  and  tabjlated  (Item  25).  If  the  site  fills,  this  is 
indicated  (item  24).  The  disposal-site  status  report  shows  for  each  period 
the  total  volume  of  consolidated  material  added  to  the  site  (item  26),  the 
resulting  volume  of  material  stored  at  the  end  of  each  period  (item  27),  and 
the  disposal  cost  (item  32).  The  elevation  (item  29)  and  surface  area  (item 

30)  corresponding  to  each  storage  are  determined  by  reference  to  the 
user-provided  capacity-elevation-area  function.  The  total  volume  of 
consolidated  material  added  and  the  cost  of  the  addition  are  presented  (item 

31)  . 


System  operation  for  the  entire  period  of  analysis  is  summarized  in  the 
dredged-material  source/disposal  site  report  (item  32).  Each 
dredged-material  source  and  disposal  site  is  shown  in  this  matrix.  The  total 
volume  moved  from  each  source  to  each  disposal  site  is  shown,  and  the  first 
and  last  periods  of  disposal  site  utilization  are  shown. 

Item  33  is  a  graphic  summary  of  disposal  site  use.  Each  period  that  any 
material  is  deposited  in  a  site  is  indicated  by  an  "X"  in  the  appropriate 
column. 
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BZAMPLB  1W0:  ANALYSIS  OP  OPIRATIMI  VITH  RISTART 


This  ezsmple  illustrates  the  "restart”  capability  of  program  D2M2.  With 
this  capability,  the  program  automatically  can  be  executed  recursively,  with 
the  final  status  of  disposal  sites  from  one  execution  serving  as  the  initial 
condition  for  the  next  execution.  The  system  considered  in  this  example  is 
the  same  as  that  considered  in  example  one,  but  instead  of  analyzing  SO 
periods  of  operation,  each  one  year  in  duration,  two  successive  2S-period 
operations  are  analyzed.  The  first  analysis  determines  the  least-costly 
system  operation  of  the  system  for  1981-2005.  Program  D2M2  then  prepares  an 
input  file  for  analysis  of  system  operation  for  2006-2030.  The  initial 
conditions  specified  in  this  input  file  are  the  final  conditions  for  2005,  as 
determined  from  the  first  analysis.  This  is  illustrated  by  Fig.  5. 

The  program  user  must  issue  appropriate  job-control  commands  to  execute 
program  D2M2  with  the  program-prepared  input  file.  The  commands  should 
accomplish  the  following:  (1)  save  the  program-prepared  restart  input  file, 
TAPEll,  from  the  first  execution;  (2)  rewind  TAPBll;  and  (3)  execute  D2I12 
with  input  read  from  TAPEll. 

Item  34  is  the  user-prepared  input  file  for  analysis  of  1981-2005 
operation.  This  input  file  is  essentially  the  same  as  the  input  file  for 
example  one.  However,  the  specified  period  of  analysis  is  25  periods  (field 
1  of  card  4).  The  value  in  field  10  of  card  4  specifies  that  D2N2  is  to 
prepare  an  input  file  to  permit  analysis  of  25  additional  periods  with 
subsequent  execution  of  the  program.  Volumes  are  specified  on  DV  cards  for 
each  dredge  site  for  only  25  periods  (cards  40-42  and  48-50).  The  volumes 
specified  for  each  site  on  DV  cards  of  the  program-generated  file  correspond 
to  the  final  value  shown  on  the  DV  cards  of  this  initial  input  file  (550,000 
volume  units  per  period  for  both  sites).  Field  3  is  blank  on  the  SX  cards 
for  the  Ocean,  Williamson  Harbor  South,  and  Williamson  Harbor  disposal  areas 
(cards  13,  19,  and  25).  This  indicates  that  these  sites  are  available 
throughout  the  original  25-period  analysis  and  the  subsequent  25-period 
analysis  performed  when  the  program  is  restarted  with  the  program-generated 
file.  The  lease  for  the  Edgewood  site  is  renegotiated  in  1999,  and  that  site 
too  is  available  throughout  the  remainder  of  the  analysis. 

Disposal-site  input  reports  and  dredge-site  input  reports  similar  to 
those  of  example  one  are  printed  but  are  not  included  here.  The  volumes 
shown  in  the  dredge  site  reports  are  those  specified  for  1981-2005. 

Item  35  is  the  capacity  expansion  iteration  log  for  analysis  of  the 
system  operation  for  1981-2005.  As  with  example  one,  this  serves  as  a 
summary  of  costs.  This  total  net  cost  shown  is  the  sum  of  acquisition,  lease 
renegotiation,  OMR,  and  operation  costs  for  1981-2005  inclusive.  The  present 
value  calculations  use  the  first  period  of  analysis,  1981,  as  the  base  period. 

System  operation  for  the  25  periods,  1981-2005,  is  summarized  by  the 
tabular  and  graphical  summaries  identified  as  items  36  and  37. 

A  restart  input  file  is  prepared  by  D2M2  to  allow  analysis  of  operation 
for  25  subsequent  periods  beginning  with  2006.  This  file  is  shown  as  item 
38.  Field  10  of  the  J1  card  (card  4)  indicates  that  1981  is  to  be  used  as 
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FIG.  5 

I 


-  Schematic  of  Restarting  Execution  of  Program  D2M2 
for  Example  Two 
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the  base  year  for  cost  amortization.  The  Initial  storage  for  each  disposal 
site  is  specified  in  field  3  of  the  SL  card  for  each  site  (cards  18,  24,  and 
30).  The  volume  to  be  dredged  at  the  Williamson  Harbor  dredging  site  and 
disposed  each  period  is  550,000  CY,  and  the  volume  to  be  dredged  and  disposed 
from  the  Williamson  H<trbor  B  site  is  550,000  (TT.  These  values,  specified 
oncards  40  and  46,  respectively,  correspond  to  the  values  previously 
specified  for  each  site  for  the  final  period  of  the  initial  analysis  (2005). 

Item  39  is  the  capacity  expansion  iteration  log  for  analysis  of  system 
operation  for  2006-2030.  The  acquisition,  renegotiation,  and  OMR  costs  shotm 
are  the  1981  present-value  equivalents  of  any  such  costs  incurred  in 
2006-2030.  The  operation  cost  shown  is  the  1981  present-value  equivalent  of 
operation  in  2006-2030.  Consequently,  the  total  net  cost  shown  is  comparable 
to  total  net  cost  computed  for  1981-2005  in  the  initial  analysis.  The 
total-net-cost  values  can  be  added  directly  to  determine  total  cost  for 
1981-2030  with  the  operation  shown. 

The  operation  for  periods  2006-2030  is  summarized  by  items  40  and  41. 

Comparison  of  example  one  results  with  example  two  results  shows  that 
the  disaggregated  analysis  may  not  yield  the  optimal  operation  scheme.  The 
cost  of  this  optimal  operation  scheme  for  50  consecutive  periods  is  $31.7 
million  while  the  sum  of  the  cost  for  the  two  consecutive  25-period 
operations  is  $32.8  million.  This  is  expected;  in  example  two,  analysis  of 
the  first  25-period  operation  cannot  consider  the  disposal  requir-/:uents  of 
the  second  25  periods.  In  example  one,  the  disposal  requirements  of  all  50 
periods  are  considered  simultaneously,  so  disposal  site  capacity  can  be 
allocated  more  efficiently.  However  the  restarting  procedure  may  provide  the 
only  means  to  analyze  long-term  operation  of  a  system  which  includes  a  large 
number  of  dredge  sites  and  disposal  sites.  Then  the  near-optimal  solution 
may  be  acceptable. 
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TOTAL  *  13750000.  *  13750000.  *  0.  *  0.  *  0.  *  0.  *  0.  *  0.  *  27500000. 
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Bzuple  Thr««:  ANALYSIS  OF  OPBBATIOM  WITH  DISPOSAL  SITB  BBSTING. 

In  practice,  dredged-materal  disposal  sites  generally  are  not  used 
continuously.  Instead  sites  are  "rested”  following  use;  the  sites  are 
allowed  to  dry  and  maintenance  is  performed  as  necessary.  This  operation 
practice  is  simulated  in  program  D2M2  with  cycles  of  availability  and 
non-availability  of  sites.  This  example  illustrates  analysis  of  system 
operation  with  site  resting.  The  system  analyzed  is  the  same  system 
considered  in  example  one,  but  the  Edgewoode  site  is  rested  for  two  periods 
of  each  five. 

Item  42  is  the  user-prepared  input  for  example  three.  This  input  file 
is  essentially  the  same  as  the  input  file  for  example  one.  The  resting  of 
the  Edgewood  disposal  area  is  specified  in  field  2  of  the  SX  card  for  that 
site  (card  31).  The  value  in  columns  9  and  10  indicates  the  site  is 
available  for  three  periods  prior  to  resting,  and  the  value  in  columns  11  and 
12  indicates  the  site  is  rested  for  two  periods.  This  cycle  is  repeated 
while  the  site  is  available  (1981-2030). 

As  with  example  one,  disposal-site  and  dredge-site  input  reports  are 
printed;  these  are  not  included  here.  Likewise,  disposal-site  status  reports 
are  printed  for  each  site;  these  are  also  omitted. 

Items  43,  44,  and  45  summarize  system  operation  for  example  three.  The 
same  volume  of  material  is  disposed  in  each  site  as  in  example  one  (see  item 
44),  but  the  timing  of  the  disposing  is  different  due  to  the  resting  of  the 
Edgewood  site  (see  item  45).  The  total  system  net  cost  is  $37.8  million 
(item  43),  which  exceeds  the  cost  of  optimal  operation  in  example  one. 
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EXAMPLE  THREE  -  D2H2  USERS  MANUAL 
ANALYSIS  OF  OPERATION  WITH  DISPOSAL  RESTING 
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Bzaapl*  Four:  AMALTSIS  OF  OPBKATION  WITH  PASTIAL  SOLUTION  SPBCIFIBD. 


Due  to  non-economlc  factors,  operation  policies  that  are  not  the 
least-costly  policies  are  sometiiaes  necessary.  These  policies  can  be 
simulated  with  program  D2H2  by  specifying  the  volume  of  material  that  is  to 
be  transported  in  a  specified  period  from  a  dredging  site  to  a  disposal 
site.  Example  four  illustrates  this. 

The  disposal  system  modeled  in  example  four  is  the  same  system  modeled 
in  example  one,  and  cards  1-S6  of  the  input  (item  46)  are  identical  to  cards 
1-56  of  the  input  for  example  1.  The  required  system  operation  is  shown  in 
Table  2.  This  operation  is  specified  on  cards  57-62.  Each  card  specifies  a 
dredged-materlal  source,  disposal  site,  period,  and  volume.  In  this  example, 
no  alternative  transporation  types  exist  between  pairs  of  dredging  sites  and 
disposal  sites,  so  transportation  types  need  not  be  specified  (field  5  of  FS 
card) .  Vfhen  a  partial  solution  is  specified,  the  user  must  specify  a 
solution  that  is  feasible:  the  disposal  site  and  dredge  site  must  be 
available  and  must  be  linked,  and  the  transportation  link  and  the  disposal 
site  capacities  must  not  be  exceeded  by  the  volume  specified. 

Disposal  site  input  reports  and  dredge  site  input  reports  similar  to 
those  of  example  one  are  printed  but  are  not  included  here.  Likewise  a 
system  volume  report  is  printed,  but  is  not  included  here. 

Item  47  is  the  capacity  expansion  iteration  log,  which  serves  here  as  a 
summary  of  system  costs.  The  present  value  of  the  operation  cost  shown  here 
represents  the  least  cost  for  operation  of  the  system,  given  the  required 
operation  specified  on  the  FS  cards. 

Items  48  and  49  are  the  components  of  the  operation  report  for  the 
Williamson  Harbor  South  disposal  site.  Similar  reports  are  printed  for  the 
other  disposal  sites,  but  these  are  not  reproduced  here.  The  volumes 
transported  to  this  site  from  the  Williamson  Harbor  dredging  site  in  1987, 
1989,  1990,  2029  and  2030  are  volumes  specified  on  the  FS  cards.  The  status 
of  the  disposal  site  with  this  operation  can  be  compared  with  the  status  of 
the  same  site  without  the  forced  operation  by  referring  to  items  22-31  of 
example  one.  In  example  one,  the  site  filled  in  2018  and  thus  was  not 
available  for  use  from  2019  to  2030.  However  in  example  four,  the 
requirement  that  550,000  be  transported  from  the  Williamson  Harbor  dredging 
site  and  disposed  at  the  Williamson  Harbor  South  site  in  both  2029  and  2030 
precludes  this.  The  disposal  site  must  have  space  available  to  store  the 
material.  This  is  accomplished  by  adding  material  to  the  site  until  2017, 
reserving  capacity  to  store  the  material  deposited  in  2019  and  2030.  The 
site  is  filled  then  in  2030. 

System  operation  for  example  four  is  summarized  in  the  dredged-material 
source/disposal  site  summary  report  and  the  disposal  site  usage  chart, 
included  as  items  50  and  51. 
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TABLE  2.  -  Required  Operation  for  Example  Pour 


Period  Dredged-material  source 

(1)  (2) 


Disposal  site 
(3) 


Volume,  in 
cubic  yards 
(4) 


1987 

Williamson  Harbor 

Williamson 

Harbor 

South 

550,000 

1989 

Williamson  Harbor 

Williamson 

Harbor 

South 

1,100,000 

1990 

Williamson  Harbor 

Williamson 

Harbor 

South 

1,100,000 

2029 

Williamson  Harbor 

Williamson 

Harbor 

South 

550,000 

2030 

Williamson  Harbor 

Williamson 

harbor 

South 

550,000 

ti^ 
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REAOIN  —  61  RECORDS  WRITTEN  TO  LOGICAL  FILE  12 


© 


u 

s 

t/i 


3 


ij} 

H* 

SjS 

s 

i/7 

=5 

i/1 

o. 

ss 

X 

UJ 

it 


iin 

,>■ 


****( 


§ 

s 


t  ♦  * 


4c  «  ^  ^  <|K  it 


<n  * 


41  4c  4C  4c  4c  4c  f 


[ 

^  c 
»—  .  ( 
o%  c 

^  C 


4c  4c  4c  4c  4c  4c  # 


4c  4c  4c  4c  4c  4c  t 

^»555 


4c  4c  4c  4c  4c  4c  4c 


iili  :i 


c 

4c  4c  4c  4c  4c  4r  c 


■  c 

flp  cp  op  flp  ‘  ( 
*  c 


«  «  4c  «  «  «  « 


71 


EXAMPLE  FOUR  -  D2H2  USERS  MANUAL 

AtMLYSIS  OF  OPEIMTION  WITH  PARTIAL  SOLUTION  SPECIFIED 
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TOTAL  *  MMOOOO.  *  3184113.  •  8126000.  *  1134474.  *  0.  •  0.  *  23066000.  *  4919247. 


DISPOSAL  SITE  STATUS  (DRY) 
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*  TOTAL  *  11533000. 


EXAMPLE  FOUR  -  02M2  USERS  MANUAL 

ANALYSIS  OF  OPERATION  WITH  PARTIAL  SOLUTION  SPECIFIED 
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Bzaapl*  riv«:  APPLICATKNI  OF  CAPACITY  EXPAMSKMi  CAPABILITY 


This  exsaple  illustrates  application  of  program  02M2  to  select  an 
economically  efficient  capacity-expansion  schedule  for  a  disposal  system. 

The  system,  illustrated  in  Fig  5,  includes  two  dredging  sites,  two  existing 
disposal  sites,  three  disposal  sites  that  will  be  acquired  in  period  2000  for 
capacity  expansion,  and  the  Williamson  Harbor  South  site,  which  may  be 
acquired  if  economically  justified.  System  operation  for  the  period 
1981-2030  is  to  be  analyzed,  but  acquisition  of  the  Williamson  Harbor  South 
site  is  considered  only  between  1981  and  2000  inclusive.  Material  from  the 
two  dredging  sites  is  removed  and  transported  to  the  disposal  sites  with  a 
16-inch  or  27-inch  pipeline  dredge  or  a  government-owned  hopper  dredge.  The 
rates  of  required  removal  from  the  dredging  site  are  the  same  as  specified  in 
example  one. 

Items  52  through  58  are  reproductions  of  selected  portions  of  the  D2M2 
output.  Item  52  is  a  list  of  the  program  input.  The  input  required  to 
exercise  the  capacity-expansion  capability  of  D2H2  is  essentially  the  same  as 
required  to  use  D2M2  for  evaluation  of  system  operation;  the  exception  is  the 
manner  in  which  the  potential  capacity  expansion  site  is  described.  As  with 
example  one,  title  cards  (cards  1-3)  and  a  J1  card  (card  4)  are  provided.  On 
this  J1  card,  the  total  number  of  disposal  sites  (field  4)  is  the  number  of 
existing  sites  plus  sites  with  specified  future  acquisition  periods  plus 
potential  capacity  expansion  sites.  A  set  of  Tl,  TD,  and  TC  cards  is 
included  for  each  dredging  and  transporting  technique.  If  techniques  not 
currently  used  may  be  used  for  expansion  sites,  these  techniques  must  be 
specified  here  with  TI,  TD  and  TC  cards.  A  set  of  SI,  SL,  SX,  SS,  SE,  and  SA 
cards  are  provided  for  each  disposal  site.  Cards  20-25  and  26-31  are  the 
disposal-site  description  cards  for  Williamson  Harbor  and  Edgewood,  the 
existing  sites.  The  three  sites  for  which  future  acquisition  is  specified 
and  the  potential  expansion  site  likewise  are  described  by  data  provided  on 
SI,  SL,  SX,  SS,  SE,  and  SA  cards  (cards  14-19,  32-37,  38-43,  and  44-49).  For 
the  capacity  expansion  site,  the  parameter  lEXP  (field  1  of  the  SX  card)  is 
three,  indicating  that  the  site  may  be  acquired  anytime  between  IPERA  and 
IPERG,  the  specified  earliest  and  latest  possible  period  for  acquisition. 

For  this  example,  these  periods  are  specified  as  1981  and  2000,  respectively. 

The  capacity-expansion  iteration  log  shows  the  expansion  plans  that  are 
selected  and  evaluated  by  the  program  and  the  results  of  each  evaluation. 

For  each  iteration,  the  trial  acqusition  period  for  the  site  is  shown  (item 
53).  For  the  potential  expansion  site,  this  acquisition  period  changes  from 
iteration  to  iteration.  The  lease  termination  period  is  identified  for  all 
sites  (item  54).  This  is  the  period  in  which  the  lease  terminates;  the  site 
then  is  no  longer  available  for  use  unless  the  lease  is  renegotiated.  For 
sites  that  are  not  leased,  this  is  shotra  as  the  period  following  the  study. 
After  the  termination  period  shown  (item  55),  the  designated  disposal  site  is 
no  longer  considered  for  use.  For  existing  or  possible  capacity-expansion 
sites  this  is  the  last  period  of  the  study.  For  leased  sites,  this  is  the 
last  period  of  the  study  if  lease  renegotiation  is  considered  in  the 
iteration;  otherwise  the  site  is  unavailable  after  the  lease  termination 
period.  Item  56  is  the  approximate  total  net  cost  of  each  alternative 
capacity-expansion  plan  evaluated.  The  estimated  present  value  of 
acquisition  cost,  lease-renegotiation  cost,  OMR  cost,  and  operation  cost  is 
included  in  this  calculation.  Acquisition  cost  is  estimated  as  a  function  of 
system  operation  in  all  cases  in  which  a  specific  acquisition  period  is  not 
defined  for  all  capacity  expansion  sites,  so  no  cost  breakdown  is  possible. 
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Williamson  Harbor 
Dredging  Site 
(WILMHA) 


Lucky  Strike  Expansion 
Site  (LUCKYS) 


Interstate  Expansion 
Site  (INTER) 


Williamson  Harbor  B 
Dredging  Site 
(WILMBO) 


Williamson  Harbor 
Disposal  Site  (WILHAR) 


Edgewood  Disposal  Site 
(EDGEWD) 


Williamson  Harbor  South 
Expansion  Site  (WILSUD) 


Belvedere  Expansion 
Site  (BELVED) 


FIG.  6.  -  Disposal  System  for  Example  Five 
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For  the  example  problem,  the  mlnimm-coat  expansion  schMw  found  calls 
for  acquisition  of  the  Williaauon  Harbor  South  site  in  1990.  The  actual 
present  value  of  total  cost  for  this  plan  is  $23786139,  which  Includes 
acquisition  cost  and  OME  cost  of  $892846  and  operation  cost  of  $22893294. 
These  true  costs  are  shown  in  the  final  row  of  the  iteration  log  (item  57), 
which  recapitulates  the  optimal  solution.  Identification  of  the  optimal  plan 
required  12  iterations  as  follows: 

Iteration  0:  Acquisition  of  the  Williamson  Harbor  South  site  is 
considered  at  any  time  following  the  earliest  specified  date  (1981). 

This  includes  not  acquiring  the  site  at  all  within  the  time  horizon  of 
this  analysis.  A  network  is  formulated  with  the  site  available  for  all 
periods  following  1981,  and  the  acquisition  cost  is  approximated  as  a 
function  of  the  utilization  of  the  site  each  period. 

Iteration  1:  At  the  beginning  of  each  iteration,  a  capacity  expansion 
site  is  selected  as  the  candidate  for  acquisition  period  revision.  In 
this  case,  only  one  potential  site  is  defined,  so  its  acquisition  period 
is  altered  from  1981-2031  to  1981-2000.  The  latest  acceptable  period 
for  acquisition  of  this  site  is  2000,  so  in  this  iteration,  the 
efficiency  of  site  acquisition  within  the  specified  allowable  time 
window  is  evaluated. 

Iteration  2:  The  utilization  of  the  Williamson  Harbor  South  site  is 
evaluated.  The  site  is  not  used  extensively,  so  acquisition  is 
postponed  until  it  is  10%  filled,  in  1997.  However  acquisition  in  lv97 
fails  to  provide  capacity  needed  earlier,  so  this  alternative  is 
infeasible . 

Iteration  3:  The  algorith’n  oack^racks  and  considers  acquisition  between 
1981  and  1996  after  determining  that  acquisition  after  1996  is 
infeasible.  The  approximate  cost  of  this  is  $23500275. 

Iteration  4:  Utilization  of  the  site  in  the  previous  iteration  is 
considered.  In  1990,  the  site  was  10%  filled,  so  acquisition  is 
postponed  until  then,  yielding  an  approximate  cost  of  $23515395. 

Iteration  5:  In  iteration  4,  the  site  never  filled  to  10%  of  capacity 
between  1990  and  1996,  so  acquisition  is  postponed  to  at  least  period 
1994.  This  is  67%  of  the  time  from  1990  to  1996,  a  percentage 
determined  from  experience.  Unfortunately,  such  postponement  is  not 
feasible  in  this  case. 

Iteration  6:  The  procedure  backtracks  and  considers  acquisition  between 
1990  and  1993.  This  is  feasible,  and  the  approximate  cost  is  $23661674. 

Iteration  7:  The  site  does  not  fill  to  10%  between  1990  and  1993,  so 
acquisition  is  postponed  using  the  67%-criteria  again.  Again  this  is 
infeasible. 

Iteration  8:  The  algorithm  backtracks  and  evaluates  the  cost  of 
acquisition  between  1990  and  1991.  The  approximate  cost  is  $23786634. 

Iteration  9:  Again  the  site  does  not  fill  to  10%,  so  again  acquisition 
is  postponed  to  1991  with  the  67%-criteria.  This  is  infeasible.  Note 
that  this  iteration  yields  a  potential  solution:  acquire  the  site  in 
1991.  However  that  solution  is  infeasible. 
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Iteration  10:  Backtracking  now  requires  evaluation  of  acquisition  in 
1990.  This  is  feasible,  and  the  cost  is  approximately  123786139.  This 
is  then  a  feasible  acquisition  plan:  acquire  the  Williaoison  Harbor  South 
site  in  1990.  However,  without  further  evaluation,  this  cannot  be 
declared  the  optimal  solution. 

Iteriation  11:  The  procedure  requires  backtracking  to  evaluate  all  as- 
yet  unconsidered  alternative  solutions.  Acquisition  between  1981  and 
1989  is  one  such  solution.  However  the  approximate  cost  exceeds  the 
cost  of  the  solution  found  in  interation  10,  so  this  solution,  and  any 
solution  requiring  acquisition  between  1981  and  1989,  can  be  eliminated 
from  further  consideration. 

Iteration  12:  The  algorithm  backtracks  finally  and  considers 
non-acquisition,  which  is  infeasible. 

Following  the  capacity  expansion  iteration  log,  the  usual  disposal  site 
reports  and  operation  summaries  (item  58  for  example)  are  printed  for  the 
system  with  the  optimal  acquisition  of  expansion  sites. 
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EXAMPLE  FIVE  -  D2N2  USERS  MANUAL 
APPLICATION  OF  CAPACITY  EXPANSION  CAPABILITY 
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APPIVDIX  IV 


D21f2  INPUT  OBSCUPTlOil 


INPUT  DESCRIPTION 


This  appendix  provides  a  detailed  description  of  program  D2M2  input 
requirements,  by  card  and  by  variable.  The  requirements  are  sumsiarized  in 
Pig.  7.  The  field  number  shoira  for  each  variable  designates  the  location  of 
that  variable  on  the  input  card.  Of  the  80  colusms  available,  coluams  1-2 
are  reserved  for  the  card  identifier  and  are  referred  to  as  field  0  (zero). 
Field  1  includes  card  colusms  3-8.  Fields  2-10  contain  eight  columns  each. 
When  necessary,  an  abbreviated  field  location  description  is  used  to  identify 
variables:  the  card  name  is  followed  by  a  decimal  point  and  the  field 
number.  For  example,  SL.8  refers  to  the  eighth  field  on  the  SL  card. 

All  variables  with  names  beginning  with  I,  J,  R,  L,  M,  or  N  are  integers 
and  must  be  right-justified  without  decimal  points.  In  general,  input 
values,  both  numerical  and  alphanumeric,  should  be  right-justif ied  in  their 
fields,  except  for  multiple  field  alphanumeric  information. 

Appropriate  values  for  each  variable  are  shown  in  the  column  headed 
VALUE.  A  *'+*'  in  the  column  indicates  a  non-negative  value  should  be 
provided.  A  **-"  indicates  a  negative  value  should  be  provided.  In  some 
cases,  either  non-negative  or  negative  values  are  acceptable;  the  variable 
description  specifies  the  interpretation  of  the  values.  The  abbreviation 
*'AN"  identifies  variables  for  which  alphanumeric  characters  are  to  be 
specified. 
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Tl,  T2,  T3 


TITLE 

These  required  cards  provide  three  lines  of  information  at  the  top  of  each 
page  of  output. 


Tl,  T2, 

T3,  CARDS 

-  TITLES 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

T1,T2, 

T3 

Card  identification. 

1-10 

TITLE 

AN 

Title  information  (center  of  title  in  card 
colusm  41) . 
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JOB  CONTSOL 


J1 


Values  on  this  required  card  define  the  period  of  analysis,  economic 
analysis  parameters,  the  dimensions  of  the  disposal  system,  and  the  output 
desired. 


J1  CARO 

-  JOB  CARD 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

J1 

Card  identification. 

1 

NPER 

+ 

Number  of  uniform  time  periods  for  which 
analysis  is  to  be  performed.  These  time  periods 
may  be  of  any  duration.  Maximum  allowable  value 
is  100. 

2 

IPERl 

+ 

Identifier  of  first  time  period  of  analysis  (for 
example,  1982).  All  benefits  and  costs  are 
discounted  to  the  beginning  of  this  time  period. 

3 

DRATE 

+ 

Effective  discount  rate  per  period,  expressed  as 
a  percentage.  This  is  used  to  determine  the 
present  value  of  all  benefits  and  costs.  A 
nonzero  value  must  be  used  because  site 
acquisition  costs  and  operation  costs  are 
amortized  in  the  analysis.  Default  value  is  IX 

4 

NDISP 

+ 

Total  number  of  existing  and  potential  disposal 
sites.  This  includes  disposal  sites  to  which 
material  is  moved  directly  from  a  dredge  site 
and  disposal  sites  to  which  material  is 
rehandled  from  other  disposal  sites.  See  note 
on  next  page. 

5 

NDREG 

+ 

Total  number  of  dredge  sites.  See  note  on  next 
page . 

6 

NTRAN 

+ 

Total  number  of  possible  dredged-material 
transportation  methods.  This  includes  dredging 
plants  (such  as  hopper  dredges)  and  rehandling 
transportation  types  (such  as  dump  trucks).  See 
note  on  next  page. 

7 

IPRNT 

0 

Transportation  link  input  report,  disposal  site 
input  report,  and  dredge  site  input  report  are 
printed  in  addition  to  list  of  input. 

1 

Transportation  link  input  report,  disposal  site 
input  report,  and  dredge  site  input  report  are 
omitted. 
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J1 


J1  CARD  (Continued) 


FIELD  VARIABLE  VALUE  DESCRIPTIOM 

8  ITNAX  +  Maximuffl  number  of  Iterations  allowed  for 

capacity-expansion  problem  solution.  Ignored  if 
system  capacity  expansion  is  not  to  be 
considered.  Default  value  is  25. 

9  NPESAM  -t-  Number  of  periods  required  for  amortization  of 

site  acquisition  or  lease  costs.  Must  not 
exceed  NPES. 

10  NPSEST  -t-  The  restart  capability  of  the  program  is  to  be 

used.  A  file  is  prepared  for  analysis  of  NPREST 
periods  following  the  last  period  of  analysis 
performed  in  this  execution  of  D2M2.  This  file 
is  a  new  input  file  for  the  same  disposal 
system.  The  initial  conditions  specified  in 
this  new  file  are  final  conditions  determined 
from  the  current  analysis.  The  file  is 
identified  as  TAPEll.  The  program  user  must 
issue  appropriate  job  control  commands  to 
re-execute  D2M2  with  this  input  file. 

0  File  for  restart  execution  of  program  is  not 

prepared. 

NOTE:  The  available  computer  memory  is  allocated  automatically  by  this 
program  based  on  the  requirements  for  storage  of  disposal  site,  dredge  site, 
and  transportation  data.  Thus  maximum  values  for  NDISP,  NDREG,  and  NTRAN  are 
not  set  explicitly.  Instead  the  following  restrictions  apply  for  the 
distributed  Harris  version  of  the  program: 

a.  (NTSAN«<34)  +  (NDISP*78)  (NDREG*18)  +  (NDREG^NPER)  +  ( (NPER+1}*12)  -f 
(NDREG’K7*MXLNKS)  +  (NPER+3)*100  must  not  exceed  40000.  MXLNRS  is  the 
maximum  number  of  disposal  sites  to  which  any  one  of  the  dredge  sites 
may  transport  material  (maximum  of  all  values  of  NTLINK  (DL.3)). 

b.  (((NDISP*2)  +  NDREG)*NPER}  +  2  must  not  exceed  5000. 

c.  (( NDISP* (2-t^N))  +  NDREG)*NPER  must  not  exceed  10000.  N  is  the  maximum 
number  of  transportation  links  that  originate  and  terminate  at  any 
disposal  site,  including  links  for  reuse  and  for  transfer. 
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TI 


TEAMSPOETATION  DKSCSIPTION 


Data  on  the  required  TI,  TD,  and  TC  cards  describe  each  possible  method 
for  moving  dredged  material  within  the  system,  including  shiiment  from  sources 
to  disposal  sites  and  shipment  between  disposal  sites  in  the  case  of 
rehandling.  NTRAN  (J1.6)  sets  of  TI,  TD,  and  TC  cards  are  required. 


TI  CARD  -  TRAMSPORTATIOM  IDBMTIFICATIOII 

The  transportation  type  identification  and  description  are  provided  on 
this  required  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KOBE 

TI 

Card  identification. 

1 

ITID 

AN 

Right-justified  6-character  transportation  type 
identifier.  This  identifier  is  used  later  to 
define  system  linkages. 

2-7 

ITNAME 

AN 

Transportation  type  description. 
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TD 


TD  CARD  -  DISTAMCBS  FOR  TRANSPORTATI(»i-COST  FUNCTION 


The  unit  cost  of  moving  material  within  the  disposal  system  is  represented 
in  02M2  as  a  function  of  distance.  The  values  on  the  required  TD  card  are  the 
distances  of  the  distance  vs.  unit  cost  function.  The  values  must  be 
specified  in  increasing  order  or  execution  will  terminate.  A  maximum  of  10 
values  may  be  provided.  Linear  interpolation  without  extrapolation  is  used, 
so  the  values  should  be  selected  accordingly.  Any  desired  units  of  distance 
may  be  used  (miles,  feet,  kilometers,  meters),  however,  the  units  of  distances 
specified  on  the  TD  card  and  the  units  of  all  other  distances  specified  must 
be  consistent. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

TD 

Card  identification. 

1 

TDIST(l) 

+ 

First  distance  value.  This  value  must  be  less 
than  the  minimum  required  transportation  linkage 
distance. 

2 

TDIST(2) 

+ 

Second  distance  value. 

3-10 

TDIST(3)- 

TDIST(IO) 

+ 

Repeat  as  above  for  values  3-10  If  less  than 

10  values  are  specified,  the  fields  remaining 
after  the  last  value  must  be  blank.  The  last 
value  specified  is  the  maximum  value  and  must 
exceed  the  length  of  any  transportation  link  in 
the  system. 
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TC 


TC  CASD  -  COSTS  FOR  TRRliSPORTATIOli-COST  FUMCTIOM 


The  values  specified  with  the  required  TC  card  are  the  unit  costs  (cost 
per  unit  volume  shipped)  corresponding  to  the  distances  specified  in  the 
corresponding  fields  of  the  TD  card.  Any  desired  units  of  cost  and  volume  may 
be  used,  but  all  costs  and  volumes  must  be  specified  with  consistent  units  of 
measurement . 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

TC 

Card  identification. 

1 

TCOST(l) 

+ 

First  unit  cost  value.  This  value  is  the  unit 
cost  associated  with  TDIST(l)  (TD.l). 

2 

TC0ST(2) 

+ 

Second  unit  cost  value.  This  value  is 
associated  with  TDIST(2)  (TD.2). 

3-10 

TC0ST(3)- 

TCOST(IO) 

+ 

Unit  cost  values  corresponding  to  TDIST(3), 
TDIST(4) . TDIST(IO)  . 
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SI 


DISPOSAL-SITB  DBSCKIPTKW  CABDS 


Data  on  these  cards  describe  each  material  disposal  site  in  turn.  NDISP 
(J1.4)  sets  of  SI,  SL,  SB,  SS,  SE,  and  SA  cards  are  required.  An  SS  card  is 
required  for  each  disposal  site  from  which  material  can  be  rehandled,  as 
identified  by  a  nonzero  value  for  NSEHND  (SB. 9). 

SI  CARD  -  DISPOSAL-SITB  IDBMTIBICATIOM 


The  dredged-material  disposal  site  identification  and  description  are 
provided  on  this  required  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

RODE 

SI 

Card  identification. 

1 

IDISP 

AN 

6-character  disposal-site  identifier. 

2-7 

IDNAME 

AN 

Disposal-site  description. 
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SL 


SL  CARD  -  DISP08AL-SZTB  LOCATION 


The  values  on  the  required  SL  card  specify  the  disposal-site  location,  the 
transportation  type  used  to  move  material  from  the  river  to  the  site,  site 
storage  data,  and  fixed  operation,  maintenance,  and  replacement  (OMR)  costs. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

RODE 

SL 

Card  identification. 

1 

SRVDST 

4 

River  distance  to  location  associated  with 
disposal  site  (see  Fig.  2).  Any  origin  for 
measurement  may  be  used,  but  all  river  distances 
must  be  specified  relative  to  the  same  origin. 

2 

SOFFMN 

4 

Distance  from  river  to  disposal  site  (see  Fig. 

2) . 

3 

IDTRAN 

AN 

Right-justified  6-character  identification  of 

the  transportation  facility  used  to  move 
material  from  the  river  to  the  disposal  site. 
This  value  must  correspond  exactly  to  one  of  the 
transportation-type  identifiers  on  TI  cards.  If 
this  field  is  blank,  IDTRAN  is  assumed  to  be  the 
same  as  the  transportation  type  used  on  the 
river  (specified  on  DT  cards),  and  SOFFMN  is 
added  to  the  river  distance  for  unit  cost 
determination.  Otherwise  the  unit  cost  is  the 
sum  of  on-river  and  off-river  unit  costs. 


4  STORl  + 


5  WDRAT  + 


6  ADDMAX  4 


0 


Volume  of  material  stored  in  disposal  site  at 
beginning  of  IPERl  (J1.2).  Must  not  exceed 
capacity  of  disposal  site. 

Average  ratio  of  volume  of  ’’wet”  material 
deposited  in  site  to  volume  of  material  after 
drying.  This  drying  is  accomplished  within  one 
time  step. 

Maximum  allowable  addition  of  unconsolidated 
material  to  disposal  site  per  period  (volume  per 
period) . 

The  absolute  value  of  ADDMAX  is  the  maximum 
elevation  increase  due  to  addition  of 
unconsolidated  material  per  time  step.  Linear 
interpolation  of  the  elevation- volume 
relationship  is  used  to  convert  this  elevation 
change  to  volume  change  for  the  network. 

The  maximum  allowable  volume  addition  per  time 
step  equals  site  capacity. 
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SL 


SL  CARD  (continued) 


FIELD  VARIABLE  VALUE 


DESCRIPTIOM 


7  STRCST  + 


8  OMRCST  + 


9  ISTPRT  0 


+ 


Cost  per  unit  volume  of  wet  material  added  to 
disposal  site.  This  cost  includes  operation 
costs  that  are  a  function  of  the  volume 
deposited  in  the  site,  including  mitigation 
costs . 

Fixed  annual  operation,  maintenance,  and 
replacement  cost  for  disposal  site.  This  cost 
is  incurred  at  the  end  of  each  period  a  site  is 
available,  regardless  of  the  volume  of  material 
stored  in  the  site. 

Disposal-site  operation  report  includes 
tabulation  of  material  added  to  and  material 
removed  from  site  each  period,  disposal-site 
status  each  period,  and  total  volumes  added  and 
removed  during  period  of  analysis. 

Disposal-site  operation  report  includes 
tabulation  of  material  added  and  removed  and 
disposal-site  status  for  period  that  site 
initially  fills,  tabulation  of  final  site 
status,  and  tabulation  of  total  voltmes  added 
and  removed. 

Disposal-site  operation  report  includes 
tabulation  of  material  added  and  removed  and 
disposal-site  status  for  period  that  site 
initially  fills  and  for  period  ISTPRT, 
tabulation  of  final  site  status,  and  tabulation 
of  total  volumes  added  and  removed. 


93 


sx 


SZ  CASD  -  nPAMSlON  SPBCIFICATION 


The  values  on  the  required  SX  card  Include  site-acquisition  data,  lease 
data,  reuse  data,  and  rehandling  data. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

SX 

Card  identification. 

1 

lEXP 

0 

This  disposal  site  is  added  to  system  at 
beginning  of  period  IPERA  (SX.3)  and  is 
available  until  lease  expires  at  end  of  period 
LLEASE  (SX.4).  The  site  is  then  deleted  from 
the  system 

1 

This  disposal  site  is  added  to  system  at 
beginning  of  period  IPERA  (SX.3).  The  lease  is 
renegotiated  at  end  of  period  LLEASE  (SX.4)  at  a 
cost  of  RENCST  (SX.5),  and  the  site  is  available 
throughout  the  analysis. 

2 

This  disposal  site  is  added  to  system  at 

beginning  of  period  IPEKA  (SX.2).  Lease 
renegotiation  at  end  of  period  LLEASE  at  a  cost 
of  SEMCST  is  considered  as  a  capacity  expansion 
option.  This  option  cannot  be  exercised  with 
the  program  "restart”  capability. 

3  Acquisition  of  this  site  is  a  capacity  expansion 

option.  The  site  may  be  acquired  between  IPERA 
and  IPESG  (SX.3),  inclusive,  at  a  cost  of  ACQCST 
(SX.IO)  and  will  then  be  available  for  the 
remainder  of  this  of  analysis.  This  option 
cannot  be  exercised  with  the  program  "restart" 
capability. 

2  IPUSE  4  Columns  9-10  contain  a  2-digit  number  defining 

the  number  of  periods  a  site  is  available  prior 
to  "resting”.  The  site  is  added  to  the  system 
in  IPERA  (SX.2)  as  specified  by  the  user  (or  as 
selected  by  the  capacity  expansion  algorithm) . 
The  site  is  then  alternately  included  for  IPUSE 
periods  and  omitted  for  IPREST  periods. 

0  Site  is  not  rested. 

IPREST  4  Columns  11-12  contain  a  2-digit  number  defining 

the  nufltber  of  periods  a  site  is  rested  in  the 
alternating  using-resting  cycles. 

0  Site  is  not  rested. 
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sx 


SX  CASD  (Continued) 


FIELD  VARIABLE  VALUE  DESCRIPTIOM 

IPERA  4  Columns  13-16  contain  a  A-digit  number.  If  lEXP 

=  0,  1,  or  2,  this  is  the  period  of  site 
acquisition.  If  lEXP  s  3.  site  acquisition  is  a 
capacity  expansion  option,  and  acquisition  will 
be  considered  only  during  or  after  IPERA.  If 
omitted,  IPERA  is  set  equal  to  IPERl  (J1.2). 

3  IPERG  4  Last  period  within  which  capacity  expansion  by 

acquisition  of  this  site  will  be  considered. 

This  value  anist  not  exceed  IPERl  4  NPER  -  NPERAM 
to  insure  that  the  site  acquisition  cost  is 
fully  amortized  within  the  period  of  analysis. 
Ignored  if  lEXP  is  not  equal  3. 

4  LLEASE  4  Period  of  lease  termination.  If  lEXP  s  0,  the 

disposal  site  is  deleted  from  the  system  at  the 
end  of  this  period.  If  lEXP  =  1,  the  lease 
renegotiation  cost  is  incurred  at  the  end  of 
this  period,  and  the  site  remains  in  the  system 
throughout  the  analysis.  If  lEXP  =  3,  LLEASE  is 
ignored.  Omit  if  site  is  available  throughout 
analysis . 

0  Lease  does  not  terminate.  Disposal  site  remains 

in  system  throughout  analysis  (including 
subsequent  "restart”  analysis,  if  that  option  is 
used) . 

5  RENCST  4  Lease  renegotiation  cost.  Required  if  lEXP  =  1 

or  2.  This  cost  is  incurred  at  the  beginning  of 

period  LLEASE4I  if  the  lease  is  renegotiated. 

6  IPERR  4  First  period  reuse  is  possible.  This  applies 

only  to  sites  with  known  acquisition  period 
(lEXP  s  0,1)  and  with  positive  value  for  REMAX 
(SX.7).  If  reuse  from  a  potential  site 
(identified  by  lEXP  =  3)  is  possible,  it  is 
assumed  to  begin  when  site  is  acquired.  In  that 
case,  this  field  is  ignored. 

7  REMAX  4  Maximum  volume  that  may  be  removed  per  period 

for  reuse.  REMAX  must  not  exceed  site  capacity. 
Omit  if  no  reuse  from  this  disposal  site. 

8  REBEN  4  Unit  benefit  of  reuse  of  material  from  this  site. 

9  NREHND  0  Material  from  this  site  cannot  be  transferred  to 

other  sites  within  the  system. 
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sx 


SZ  CAID  (Continued) 


riBLD  VARIABLE  VALUE  DESCRIPTIOM 

4-  Number  of  sites  to  which  sisterial  from  this  site 

may  be  transferred.  Maximum  value  is  2.  Data 
describing  the  transportation  links  to  these 
sites  are  provided  on  the  SR  card. 

10  AC()CST  -f  Total  site  acquisition  cost.  This  cost  may 

include  cost  of  facilities,  of  real  estate,  and 
of  mitigation  of  adverse  environmental  or 
cultural  impacts. 
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SR  CARD  -  TRANSFER  DESCRIPTION 


SR 


The  data  on  this  optional  card  describe  the  potential  facilities  for 
transferrins  material  from  the  disposal  site  identified  on  the  SI  card.  The 


SR  card 

is  required 

if  NREHND 

(SX.9)  is  positive. 

FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

SR 

Card  identification. 

1 

IRESIT(l) 

AN 

Right-justified  6-character  identification  of 
the  first  disposal  site  to  which  material  may  be 
transferred.  This  site  must  be  described  also 
by  data  on  SI,  SL,  SX,  SR,  SS,  SE,  and  SA  cards. 

2 

IREPER(l) 

+ 

First  period  that  transferring  from  site 
IDISPiSI.l)  to  site  IRESIT  (SR.l)  is  possible. 
This  applies  only  to  sites  with  known 
acquisition  period.  Otherwise  transfer  is 
assumed  to  begin  when  both  sites  IDISP  and 

IRESIT  have  been  acquired. 

3 

REDISK  1) 

+ 

Distance  to  site  IRESIT. 

4 

IRETYPd) 

AN 

Right-justified  6-character  identification  of 
the  transportation  method  for  transferring 
material.  This  value  must  correspond  to  one  of 
the  transportation  types  specified  previously  on 
TI  cards . 

5 

REHHAX(l) 

+ 

Haximum  volume  per  period  that  may  be 
transferred  to  site  IRESIT.  If  no  value  is 
specified,  maximum  site  capacity  is  used. 

6 

IRESIT(2) 

AN 

Same  as  IRESIKl),  but  for  second  transfer  site, 
if  one  is  available. 

7 

IREPER(2) 

+ 

Same  as  IREPERCl),  but  for  second  transfer  site, 
if  one  is  available. 

8 

REDISK  2) 

+ 

Same  as  REDISKl),  but  for  second  transfer  site, 
if  one  is  available. 

9 

IRETYP(2) 

AN 

Same  as  IRETYPd),  but  for  second  transfer  site, 
if  one  is  available. 

10 

REHMAX(2) 

+ 

Same  as  E^HMAXd),  but  for  second  transfer  site, 
if  one  is  available. 
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SS,  SB 


SS  CARD  -  DISPOSAL-SITB  CAPACITY 


The  data  on  the  required  SS,  SB,  and  SA  cards  describe  the  capacity, 
surface  elevation,  surface  area  relationship  of  the  disposal  site.  Any 
consistent  system  of  measurement  may  be  used.  The  values  are  specified  such 
that  CAPCTY(l)  is  the  capacity  corresponding  to  ELEV(l),  «rhich  corresponds  to 
AREA(l).  The  storage  values  specified  must  increase  or  execution  will 
terminate.  The  same  number  of  values  must  be  specified  on  each  of  the  cards. 
A  maximum  of  9  values  may  be  provided.  Linear  interpolation  without 
extrapolation  is  used,  so  the  values  specified  should  be  selected 
accordingly.  The  last  nonzero  value  of  the  SS  card  is  assumed  to  be  the 
maximum  allowable  volume  of  aiaterial  deposited  in  the  site. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

SS 

Card  identification. 

1 

Not  used. 

2 

CAPCTY(l) 

+ 

First  storage  value 

3-10 

CAPCTY(2)- 

CAPCTY(9) 

+ 

Storage  values  2-9. 

SB  CARD  -  DISPOSAL-SITE  ELEVATION 


The  required  SE  card  identifies  the  elevations  corresponding  to  the 
capacity  values  of  the  SS  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

SE 

Card  identification. 

1 

Not  used. 

2 

ELEV(l) 

+ 

First  elevation  value.  This  is  the  elevation 
that  corresponds  to  CAPCTY(l). 

3-10 

ELEV(2)- 

ELEV(9) 

+ 

Elevation  values  2-9,  corresponding  to  storage 
values  2-9  on  SS  card. 
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SA 


SA  CARP  -  DISPOSAL-SITE  AREA 


The  surface  area  values  of  the  capacity-elevation-area  relationship  are 
specified  on  this  required  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KOBE 

SA 

Card  identification. 

1 

Not  used. 

2 

AREAvl) 

+ 

First  surface  area  value.  This  is  the  surface 
area  that  corresponds  to  CAPCTY(l). 

3-10 

AREA(2)- 

AREA(9} 

+ 

Area  values  2-9,  corresponding  to  storage  values 
2-9  on  SS  card. 
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DI 


OEKDGB-SITB  DBSCBIPTION  CARDS 


Data  on  these  required  cards  describe  each  dredge  site  in  turn.  NDREG 
(J1.5)  sets  of  DI,  DL,  DT,  and  DV  cards  must  be  included.  Within  each  set, 
one  DT  card  is  included  for  each  possible  linkage  of  the  dredge  site  with  any 
disposal  site. 


DI  CARD  -  DRBD6B-SITB  IDBMTIFICATION 


The  dredge-site  identification  and  description  are  specified  on  this 
required  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

RODS 

DI 

Card  identification. 

1 

IDREG 

AN 

Right-justified  6-character  dredge-site 
identifier. 

2-7 

IDRNAM 

AN 

Dredge  site  description. 
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DL 


DL  CASD  -  DSBDGB-SITE  LOCATION 


The  data  on  the  required  DL  card  specify  the  dredged-material  source 
location,  the  nunJjer  of  links  of  the  source  with  disposal  sites,  and  the 
multiplier  that  may  be  used  to  scale  all  the  specified  material  volumes  for 
this  source. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

DL 

Card  identification. 

1 

DRVDST 

+ 

River  distance  to  dredged-material  source 
location.  The  origin  of  measurement  must  be 
constant  and  the  same  used  for  specification  of 
disposal-site  location. 

2 

DOFFMN 

+ 

Distance  from  river  to  dredge  site.  If  dredge 
site  is  on  a  tributary,  this  is  the  distance 
from  the  point  of  intersection  of  the  tributary 
and  the  main  channel. 

3 

NTLINK 

+ 

Number  of  possible  links  of  this  dredged 
material  source  with  existing  or  potential 
disposal  sites.  These  links  are  defined  by  DT 
cards  that  follow. 

4 

VMULT 

+ 

Material- volume  multiplier.  All  material 

volumes  specified  on  DV  cards  for  this  dredge 
site  are  multiplied  by  VMULT  prior  to  analysis. 


0  VMULT  is  assumed  equal  to  1.0. 
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DT 


DT  CARD  -  DRSDGE  SITE  TO  DISPOSAL  SITE  TEAESPORTATION  DESCRIPTION 


One  DT  card  is  required  to  describe  each  possible  link  of  this  dredged- 
material  source  to  existing  or  potential  disposal  sites.  A  total  of  NTLIMK 
(DL.3)  cards  are  required. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

DT 

Card  identification. 

1 

LDID 

AN 

Right-justified  6-character  identifier  of  site 
to  which  material  may  be  transported  for 
disposal.  This  identifier  must  match  exactly 
one  previously  specified  on  SI  cards. 

2 

LMTYP 

AN 

Right-justified  6-character  identifier  of  method 
used  for  transportation  of  material  on  river. 
This  identifier  must  match  exactly  one 
previously  specified  on  TI  cards. 

3 

LOFTYP 

AN 

Right-justified  6-character  identifier  of  method 

used  to  transport  dredged  material  from  source 
to  river.  This  identifier  must  match  exactly 
one  previously  specified  on  TI  cards.  Omit  if 
DOFFMN  (DL.2)  is  zero.  If  this  field  is  blank, 
LOFTYP  is  assumed  to  be  the  same  as  LMTYP 
(DT.2).  In  that  case,  DOFFMN  (0L.2)  is  added  to 
the  river  distance  for  unit  cost  determination. 

4  VLHAX  +  Maximum  allowable  volume  per  period  that  may  be 

transported  via  this  transportation  link.  If  no 
value  is  specified,  the  entire  volume  removed 
from  this  site  each  period  is  used  as  the  limit. 
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DV 


DV  CAKD  -  DBIDCID-HATKRIAL  VOLUHKS 


The  volume  to  be  transported  from  the  dredge  site  each  period,  before 
drying,  is  specified  with  the  required  DV  card.  Ten  values  are  specified  per 
card.  NPER  (Jl.l)  values  must  be  specified,  unless  a  constant  volume  is  to  be 
used.  In  that  case  the  value  specified  for  the  first  period  is  used  in  all 
subsequent  periods.  This  is  controlled  by  the  presence  of  -1  in  field  3.  If 
VMULT  (DL.4)  is  nonzero,  all  volumes  specified  or  implied  are  multiplied  by 
VMULT. 


FIELD  VARIABLE  VALUE 


DESCRIPTIOil 


0  RODE  DV 


Card  identification. 


1 


Not  used. 


2  VOLMATd)  4 


Volume  of  ”wet**  material  to  be  transported  from 
this  dredge  site,  period  IPERl. 


3  VOLMAT(2)  + 


Volume  of  ”wet**  material  to  be  disposed,  period 
IPERl+1. 


-1 


4  VOLMAT(3)  + 


5-  VOLMAT(4)-  + 

NPER  VOLMAT 

(NPER) 


The  volume  of  material  for  period  IPERl^l  and 
for  all  subsequent  periods  is  VOLMATd) ^^VMULT. 
No  additional  values  are  required  on  this  card 
and  additional  DV  cards  for  this  site  are  not 
required. 

Volume  of  ”wet”  material  to  be  disposed,  period 
IPER142.  Omit  this  value  and  subsequent  values 
if  VOLMAT(2)  =  -1. 

Volume  of  ”wet”  material  to  be  disposed  in 
remaining  periods.  If  NPER  exceeds  9,  use 
addition  cards  as  needed,  with  9  values  per 
card,  beginning  with  Field  2  of  each  card. 


103 


FS 


USBR-SPBCIFIEO  SOLUTION 


Data  on  these  optional  cards  force  transportation  of  a  specified  volume  of 
dredged-material  from  a  specified  source  to  a  specified  disposal  site 
regardless  of  the  cost.  If  multiple  transportation  linkages  exist,  the 
transportation  type  must  also  be  specified.  The  user  must  realize  that  fixing 
a  portion  of  the  solution  in  this  manner  may  yield  a  network-flow  programning 
problem  to  which  no  solution  exists,  especially  in  the  case  of  capacity 
expansion  sites.  Furthermore,  the  user  must  define  a  feasible  solution;  the 
linkage  must  exist,  the  volume  to  be  transported  must  not  exceed  the  capacity 
of  the  disposal  site,  and  transporting  and  disposing  the  material  must  not 
cause  overfilling  of  the  disposal  site  in  a  single  period. 


FS  CARD  -  FIXED  SOLUTION  SPECIFICATION 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

KODE 

FS 

Card  identification. 

1 

JDREG 

AN 

Right-justified  6-character  identifier  of  the 
dredged-material  source.  This  value  must 
correspond  exactly  to  one  of  the  sources 
identified  on  DI  cards. 

2 

JDISP 

AN 

Right-justified  6-character  identifier  of  the 
disposal  site  to  which  material  is  to  be 
transported.  This  value  must  correspond  exactly 
to  one  of  the  disposal  sites  identified  on  SI 
cards . 

3 

JPR 

+ 

Period  in  which  specified  value  is  to  be 
transported.  Both  IDREG  and  IDISP  must  be 
available  in  period  IPER  or  execution  will 
terminate. 

4 

VOL 

+ 

Volume  to  be  transported.  Must  not  exceed 
volume  specified  on  DV  card  for  source  IDREG, 
period  IPER. 

5 

JTRTYP 

AN 

Right-justified  6-character  identifier  of  the 
transportation  method  used,  if  more  than  one 
type  of  linkage  exists.  This  value  must 
correspond  exactly  to  one  of  the  transportation 
types  identified  on  DT  cards. 
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FIG.  7.  -  Summary  of  Program  D2M2  Input 
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